MINIMALITY OF THE DATA IN WAVELET FILTERS 
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with an Appendix by Brian Treadway 



Abstract. Orthogonal wavelets, or wavelet frames, for (R) are associated 
with quadrature mirror filters (QMF), a set of complex numbers which relate 
the dyadic scaling of functions on R to the Z-translates. In this paper, we 
show that generically, the data in the QMF-systems of wavelets is minimal, in 
the sense that it cannot be nontrivially reduced. The minimality property is 
given a geometric formulation in the Hilbert space (Z) , and it is then shown 
that minimality corresponds to irrcducibility of a wavelet representation of the 
algebra O2; and so our result is that this family of representations of O2 on the 
Hilbert space (Z) is irreducible for a generic set of values of the parameters 
which label the wavelet representations. 



1. Introduction 

Let L"^ (R) be the Hilbert space of all L^-functions. For G L'^ (R), set 

(1.1) %l)n,k {x) 2^ V (2"a; - k) for x E R, and n,ke'Z. 

We say that ^p is a wavelet (in the strict sense) if {'ipn,k ; n, k € 1j} constitutes an 
orthonormal basis in (R); and we say that ^ is a wavelet in the frame sense 
(tight frame) if 

(1-2) uwUm)- E i(v^",^i/>i' 

holds for all / £ (R), where (• | •) is the usual (R)-inncr product, i.e., 
{4'n,k I /) = Jb 'ipn.k jx) f jx ) dx — Cn,k- The uumbers c„^k are the wavelet coeffi- 
cients. It is known |Dau92 , H6r95| that a given wavelet ip in the sense of frames is 



a (strict) wavelet if and only if ||V'|Il2(r') = 1. We shall have occasion to consider 
scaling on R other than the dyadic one, say x 1-^ Nx where N £ N, N > 2. Then 
the analogue of (^]^) is 

(1.3) 'ipn,k{x) --^ N^i(;{N"-x - k) , x£R,n,k£Z. 



However, in that case, it is generally not enough to consider only one ^|J in L 



2 , 



If the wavelet is derived from an iV-subband wavelet filter as in |BrJoOC|; then we 
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construct ^p''^\'^p^'^\ . . . ,"0'^ in such that the functions in ( |l.3| ) have the 

basis property, either in the strict sense, or in the sense of frames. Then the system 



(1.4) 



l<i<N, 



n,k € Z I 



constitutes an orthonormal basis of (M), or, alternatively, a tight frame, as in 
( |l.2| ) but with the V'n i functions in place of ijjn.k- 

Our main point is to show how the notion of irreducibility for representations of 
the Cuntz algebra On corresponds to optimality of the corresponding wavelet filters. 
Since we are addressing two different audiences (wavelets vs. representation theory), 
a few more details are included in this paper than might otherwise be customary. 
Our main result is that the irreducibility of the representation (equivalently, mini- 
mality of the filter) is generic for the wavelet representations; see Theorems 5.9 and 



3.7. In addition, we show that generically, two different filters yield inequivalent 
representations, i.e., the corresponding two representations are not unitarily equiv- 
alent. This was known earlier only in very restrictive special cases [BrJoOOj, and 
the general case treated here has not previously been discussed in the literature. 
Moreover, the methods used for the special cases in fact do not at all carry over 
to the general case. We are concerned with the wavelet filters which enter into 
the construction of -f/;'-'^' , -f/''^'' 



,1/''^"^^ in (|Tj). These filters (see (|l|)-(|L3) and 
( ^/7| ) below) are just a finite set of numbers which relate the Z-translates of these 
functions to the corresponding scalings by a; i-^ Nx. Hence the analysis may be 
discretizcd via the filters, but the question arises whether or not the data which go 
into the wavelet filters are minimal. Representation theory is ideally suited to make 
the minimality question mathematically precise. (This is a QMF-multiresolution 
construction, an d it is it s minimal it y and e ffi ciency wh ich c oncern us here. While it 
is true, see, e.g., iGab98| , [fPT9|, lBagOO| , [|BaMe99t| , and ||DaiL98| , that there are 



other and different possible wavelet constructions, it is not yet clear how our present 
techniques might adapt to the alternative constructions, although the approach in 
|DaiL98 | is also based on operator-theoretic considerations.) 

To explain the minimality issue for multiresolution quadrature mirror (QMF) 
wavelet filters, we recall the scaling function </? of a resolution in (R). Let g G N, 
and let oq, oi, . . . , a2g-i be given complex numbers such that 



(1.5) 

and 
(1.6) 



2g-l 
k=0 



E 



o.k+2iak — 



if / = 0, 
■lil^O. 



In the summation (^|^), and elsewhere, we adopt the convention that terms are 



L6D 



defined to be zero when the index is not in the specified range. Hence, in 
it is understood that ak+21 = whenever k and I are such that k + 21 is not 
in {0, l,...,2.g-l}. It is known |BrJoOO[ |BEJO0|, |Mal99| that there is a c/j e 



(M) \ {0} of compact support, unique up to a constant multiple, such that 



(1.7) 



2g-l 

ip{x) ^ E "'"'^ " 
k=0 



k), 
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in fact, supp {(p) C [0, 2g ~ 1]. (If H denotes the Hilbert transform of (R), and (p 
solves (^), then Hip does as well; but H(p will not be of compact support if ip is.) 
In finding ip in (1.7), there are methods based on iteration (see Appendix), on ran- 
dom matrix products, and on Fourier transform, see | BrJoOC| , [BEJOOj, | BrJo99bt , 
|Coh92b| , iCoRy95| , and pau92| ; and the various methods intertwine in the anal- 
ysis of ip, i.e., in deciding when ip (x) is continuous, or not, or if it is diffcrcntiable. 
This issue will be resumed in the Appendix below, which is based on |BrJo99b|. 
But the next two sections will deal with the minimality question alluded to above. 

1.7), and let Vq be the closed subspace in Ti (:= (R)) spanned 
, by the integral translates of the scaling function ip. Let 



Let ip be as in 
by {(/? (x - fc) ; fc eZ} 
U {:= Un) be 



(1.8) 



Uf (x) ■.= N~2f 



N 



the unitary scaling operator in H = (M). Then if TV = 2, 

(1.9) UVo C Vo 
is a proper subspace, and 

(1.10) /\c/«Vo = {0}; 



see fiBEJOOl and [|Dau92| , Ch. 5]. Setting V„ := f7"Vo and 
(1.11) W„:=V„_ieV„, 



we arrive at the resolution 
(1.12) 



n>l 

and the wavelet function ^ is picked in Wq ; see Table |l]. We will set up an iso- 
morphism between the resolution subspace Vo and l'^ (Z), and associate operators 
in (Z) with the wavelet operations in Vo C (M) . This is of practical signifi- 
cance given that the operators in (Z) are those which are defined directly from 
the wavelet filters, and it is th e digital filter operations which lend themselves to 
algorithms. Generalizing ( 1.11 ), in the case of scale N (> 2) the space Vo UnVq 

1; see (|3l^)-(|37l7|). 



splits up as a sum of orthogonal spaces yV{ 



(0 



1,2,. 



2. Representations of On and Table |i] (Discrete vs. continuous 

wavelets) 

The computational significance of the operator system in Table |l| (scale N — 2) 
is that the operators which generate wavelets in (M) become modeled by an as- 
sociated system of operators in the sequence space i"^ (:= P (Z) ^ (T)). (We will 
do the discussion here in Section ^ just for N ^ 2, but this is merely for simplicity; 
it easily generalizes to arbitrary N.) Then the algorithms are implemented in i'^ by 
basic discrete operations, and only in the end are the results then "translated" back 
to the space L^(R). The space L^(R) is not amenable (in its own right) to discrete 
computations. This is made precise by the frame operator W : (= (T)) Vo 
(C (R)) defined as 



(2.1) 



W:e^3 (6) ^ i^-k)^ 



feez 
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Table 1. Discrete vs. continuous wavelets, i.e., i"^ vs. (M) 



{0} 



V2 



Vi 



Vo 



\ 



\ 



finer scales 



W2 



Wi 



rest of L 



W 



{0} 



u 


u 


U 


So 


1 

So 


1 

So 



w 



Sic 



SoJC 



C = SiP 



2p2 



S^i 



So^ 



If has orthogonal translates, then W will be an isometry of onto Vo (c (R))- 
Even if the functions {ip {x — fc)}^^^ formed from Lp by Z-translates only constitute 
a frame in Vq, then we will have the following estimates: 



(2.2) 



c\''-m\,.<\\w£,\\ 



< c. 



1/2 



leii. 



where ci and C2 are positive constants depending only on ip. 

Lemma 2.1. // the coefficients {ofc ; fc = 0, 1, . . . , 2^ — 1} from (1.7) satisfy the 
conditions in ( |l.6| ) , then the corresponding operator Sq : 

(2.3) (S'oOfe = ^ ^ ^ a/c-2/'?; = ^ >^ flp^fc-p, 



pGZ: 
p=k mod 2 



£ , given by 

fc e z. 



is isometric otic? satisfies the following intertwining identity: 

(2.4) WSo = UW, 

where U is the dyadic scaling operator in (R) introduced in (1.8). (Here we 
restrict attention to N — 2, but just for notational simplicity!) Setting b^ :— 
(— l)'^ a2g-i-fe, and defining a second isometric operator Si: l"^ by formula 

( ^.3| ) with the only modification that (bk) is used in place of [ok), we get 

(2.5) S*S^. = (5j,fcl£2 
and 

(2.6) 



which are the Cuntz identities from operator theory |Cun77|, and the operators So 
and Si satisfy the identities indicated in Table 0. 

Remark 2.2. For understanding the second line in Table ||, note that S'o is a shift 
as an isometry, in the sense of [ 3zFo7C | , and C := Sii"^ is a wandering subspace for 
So, in the sense that the spaces £, Sq£, SqC, . . . are mutually orthogonal in i"^. 
To see this, note that (U) imphes that (£ :=) Sii"^ =fQ Soi^ = ker (S^). As a 
result, we get the following: 
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Corollary 2.3. The projections onto the orthogonal subspaces in the second line 
of Table 1 corresponding to the Wi, W2, . . • subspaces of the first line {see (1.11)) 



proj£ = SiSl = I - SqSq, 



proj s^-'c = 5r '5*"-^ - sssr 



Proof. Immediate from Lemma p.l|, Remark 2.2, and (1.12) 



□ 



Remark 2.4. Any system of operators {Sj} satisfying (2.5)-(2.(;) is said to be a 
representation of the C*-algebra O2, and there is a similar notion for On when 
N > 2, with On having generators Sq, Si, ... , Sn-i, but otherwise also satisfying 
the operator identities (2.5)-(2.(: ). The power and the usefulness of the multires- 



olution subband filters for the analy sis of wavelets and their alg orithms w as first 
demo nstrate d forcefully in CoWi93 and Wic93 |; se e es pecially | CoWi93 , p. 140] 
and I Wic93| , p. 157], where the Cjy-relations (^.5|)-(2.6) a re ident ified, and ana- 
lyzed in the case N = 2. Around the same time, A. Cohen |]Coh92b l identified and 
utilized the interplay between £^ and (M) which, as noted in Section || above, is 
implied by the OAr-relations and their representations. But neither of those prior 
references takes up the construction of OAr-representations in a systematic fashion. 
Of course the quadrature mirror filters (QMF's) have a long history in electrical 
engineering (speech coding problems), going back to long before they were used 
in wavelets, but the form in which we shall use them here is well articulated, for 
example, in []CEG77 . Some more of the history of and literature on wavelet filters 
is covered well in [Mey93| and | BenO0[ . 



Definition 2.5. A representation of On on the Hilbert space i'^ is said to be irre- 
ducible if there are no closed subspaces {0} ^ T i p $ l'^ which reduce the represen- 
tation, i.e., which yield a representation of (2^)-(2_^) on each of the two subspaces 
in the decomposition 



(2.7) 



f ^Ho® {f e Ho 



where e Ho = (Hq)^ = {£. e £^ : = 0, V?/ G Hq}. 



Pro of of Lem ma |2.1| . Most of the details of the proof are contained in | BrJo97b) 
and [ BrJoOC ], so we only sketch points not alr ead y covered there. The essential step 
(for the present applications) is the formula (2^), which shows that W intertwines 
the isometry S'o with the restriction of the unitary operator U : f ^ -^f {x/2) to 
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the resolution subspace Vo C (M). We have: 
(UWO (i) = -L(H'f)(|) 

= 7f XI XI ^'^"''^ (a; - 2fc - 
- (M^SoO {x) 



(by (ii|)) 



(by (0)) 



(by (P)) 



(by dgj)) 



for aU ^ e j^ii X e M. This proves (2.4). 



□ 



For later use, we record the operators on the respective Hilbert spaces (T) = i"^ 
and (R), and the corresponding transformation rules with respect to the operator 



W . Let N be the scale number, and let {ak)klo ^ be given satisfying 



(2.8) 



and set mo (z) S*^o ^ ^kZ^ ^ z ^T. The following summary table of trans- 

formation rules may clarify the proof. 
(2.9) 

TRANSLATION 



SCALING 

l2(R): f^-^f(^) F{x)^F{x~l) 
t w 

: f ^ 5]afe-Ar/6 (a) ^ (a-i) 

Fourier transform ' 

L2(t): /^^o(^)/(^^) /(^)^^/(^) 



real wavelets 
discrete model 



periodic model, 

T = M/27rZ 



Remark 2.6. The significance of irreducibility (when satisfied) is that the wavelet 
subbands which are indicated in Table |l| are then the only suhhands of the corre- 
sponding multiresolution. We will show that in fact irreducibility holds generically, 
but it does not hold, for example, for the Haar wavelets. In the simplest case, the 
Haar wavelet has g = 2 = N , and the numbers from Lemma 2.1 are 

(2.10) "'^ ^ 



ao 
bo 



ai 
bi 



1 



1 



Hence, for this representation of O2 on we may take TCq — i"^ (0, 1,2,...), and 
therefore Hq = £^ (. . . , — 3, — 2, — 1). Returning to the multiresolution diagram 
m Table |l|, this means that we get additional subspaces of (R), on top of the 
standard ones which are listed in Table |l|. Specifically, in addition to 

V„ = [/"Vo = WS^'f and >V„ = V„_i 9 V„ = WSJ^-^S^f, 
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we get a new system with "twice as many", as follows: V,!^'' and yvi^\ where 
V(+) = WSJ} (Ho) , = WSJ^-'S, {Ho) ; 



and 



For the case of the Haar wavelet, see ( 2.10 ), 

V,$+^ C L^(0,c3o), Vt'> C L'^ (-00,0) , 
or rather, Vo consists of finite linear combinations of Z-translates of 



(2.11) 



(p{x) 



1 if < X < 1, 
if a; e R\ [0,1), 



i.e., functions in (R) which are constant between n and n + 1 for all n G Z; and 
(2.12) V^+^ ==Vo 2.2(0, cx)), V^"^ Vo n (-00, 0) . 

Hence we get two separate wavelets, but with translations built on {0, 1,2,...} and 
{. . . , —3, —2, —1}. In view of the graphics in the Appendix below, it is perhaps 
surprising that other wavelets (different from the Haar wavelets) do not have the 
corresponding additional "positive vs. negative" splitting into subbands within the 
Hilbert space (R). 



Remark 2.7. There are other dyadic Haar wavelets (mock Haar wavelets), in ad- 
dition to (2.11). For example, let 



(2.13) 



\/2fc+l 







if < a: < 2fc + 1, 
if a; e R \ [0, 2fc + 1) 



Then it follow s tha t there is a splitting of Vo into orthogonal subspaces which is 
analogous to ( 2.12 ), but it has many more subbands than the t wo, "positive vs. 
negative", which are special to the standard Haar wav elet (^.11 ). For details on 
these other Haar wavelets, and their decompositions, see [BrJo99a , Proposition 8.2]. 
They are only tight frames, and the m-functions of ( |2.13|) are 



(2.14) mo (z) = ^ (1 + ^''^') : ^1 (^) = 7f (1 - 



zeT. 



Hence, after adjusting the C'2-representation T with a rotation V G \J2 (C), we have 
(2.15) To/(z) = /(z2), T,f{z) = z''^+^f{z'), f e L^T) ^ £^ 



and the two new operators Tq,Ti will satisfy the C'2-identities (2.5)-( |2.6| ); the 
representation will have the same reducing subspaces as the one defined directly 
from niQ and mi. The explicit decomposition of the multiresolution subspaces 
corresponding to (2.12) may be derived, via W in Table 0, from the decomposition 
into sums of irreducibles for the C'2-representation on P which corresponds to 
( 2.12 ). This means that the decomposition (2/7) associated with ( 2.15 ) and ( 2.15 ) 
has more than two terms in its subspace configuration. 
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3. Wavelet filters and subbands 

The operators of wavelet filters may be realized on either one of the two Hilbert 
spaces £^(Z) or L^(T), T = R/27rZ, and L^(T) defined from the normalized Haar 
measure fj, on T. But, of course, £^(Z) = L^(T) via the Fourier series. For a given 
sequence oq, ai, . . . , UNg-i, consider the operator Sq in £^(Z) given by 



(3.1) 



^ I — > So^ and {So^k 



Setting mo(z) = -j= Y,k=o ^ "-^^^ ^^'^ 



(3.2) 



^ ^k=0 O'kZ 



/eL2(T), zeT, 



we note that 5*0 and Sq are really two versions of the same operator, i.e., that 
{Sof)^= Soif) when / — (^fc) from the Fourier series. (The first one is the discrete 
model, and the second, the periodic model, referring to the diagram (^.9|).) Hence, 
we shall simply use the same notation Sq in referring to this operator in either one 
of its incarnations. It is the (3.1) version which is used in algorithms, of course. 
Let ip £ i^(R) be the compactly supported scaling function solving 



(3.3) 



fix) 



Ng-l 

akip{Nx — k). 

k=0 



Then define the operator W : e^{Z) L^{m.) by (^. The conditions on the 
wa velet filter {ak} in (L5)-(L6) and (2^) may now be restated in terms of 7719(2:) 
in (3.2) as follows: 



(3.4) 

and 
(3.5) 



AT-l 

E 

k=0 



mo{ze'- «") 



mo(l) = VN. 



It then follows from Lemma 2.1 that W in ( ^.l| ) maps ^^(Z) onto the resolution 
subspace Vo (c L^(R)), and that 



(3.6) 



UnW = WSo 



where UNf{x) = N'^'^/'^ f {x/N), f e i^(M), x e R. We showed in |BrJoOO| that 
there are L°°-functions mi, . . . ,tojv-i such that the N-hy-N complex matrix 



(3.7) 

is unitary for all z G T. If we define 
(3.8) ^,/(z)-m,(z)/(z^) 
then 



N-l 



j,k=0 



(3.9) 



SjSk — 5j,klL^{Y), 
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and 

(3.10) ^5,5; = /i.(T). 

((|3.5|) is not needed for this, only for the algorithmic operations of the Appendix. 



Lemma 3.1. The solutions (™j)j^o ^'^ (3-7) are in 1-1 correspondence with the 
semigroup of all polynomial functions 

(3.11) A:T — ^UAr(C), 

where \]pf{C) denotes the unitary N x N matrices. 

Proof. The correspondence is to A with 



N-l 



(3.12) 



,(z)=^A,,fc(z^). 



fc=0 



□ 



and in the reverse direction, 

(3.13) A,^k{z) = ^ E ^^'rn.iw) 

does the job, as can be checked by direct substitution. 

We also showed in [ BrJoOOf that if toq is given, and if it satisfies ( |3.4[ ), then it is 
possible to construct toi, . . . , nriN-i such that the extended system toq, toi, . . . , to^t-i 
will satisfy (3.7). As a consequence, A in ( |3.13 ) will be a UAr(C)-loop, and the orig- 
inal Too is then recovered from (3.12) for j = 0. To stress the dependence of the 
operators in ( |3.8| ) on the loop group element A we will denote the corresponding 
operators T^^\ and it follows that, \i A= l^r, then the operators Si of (3.i 

z*/ (z^) , where i = 0, 1, 



are 



(3.14) 



,A^- 1, 



and we will reserve the notation Si for those special ones, i.e.. Si 

Let Sj ^ t'^-^^ be an arbitrary wavelet representation. By virtue of (3.£)-(3.10), 
i2(T), or equivalently splits up as an orthogonal sum 

(3.15) T\^\e{Z)), J = 0,1,..., TV -1. 



We saw that the wavelet transform W of (2J) maps ^■^(Z) onto Vo, and from ( |3.6|) 
we conclude that W maps T^'^\P{Z)) onto Un{Vo) {-. Vi). Hence, in the A^- 
scale wavelet case, W transforms the spaces T\^\f-{'L)) (C ^2(Z)) onto orthogonal 



subspaces 

(3.16) 

where 
(3.17) 



0) 



, Af- 1 in L'^i 



3 

and 



;0) 



J = 1, . 



,Af-l. 



Each of the spaces Vi and Wl is split further into orthogonal subspaces corre- 
sponding to iteration of the operators Tg'^\T^'^\ . . . , T^'^ of (3.9)-( ^J^ ). It is the 
system {Tj^^}^Sq^ which is called a wavelet representation, and it follows that the 



10 



PALLE E. T. JORGENSEN 



wavelet decomposition may be recovered from the representation. Moreover, the 
variety of all wavelet representations is in 1-1 correspondence with the semigroup of 



polynomial functions A in ( |3.1l[) . Operators {Tj '} satisfying (|3.9[)-( |3.10| ) are said 



to constitute a representation of the C*-algebra On, the Cuntz algebra [Cun77|, 



and it is the irreducibility of the representations from (p.q) which will concern us. 



If a representation (3.8) is reducible (Definition 2.5), then there is a subspace 
(3.18) 0CHoCl2(t) 

which is invariant under all the operators Tj^^ and T^j"^^ * , and so the data going 
into the wavelet filter system {rrij} are then not minimal. 

4. A LEMMA ABOUT PROJECTIONS 

Our main result is that for a generic set within the class of all wavelet represen- 
tations, we do have irreducibility, i.e., there is no reduction as indicated in ( |3.18 ) in 



Section ^ In proving this, we will first reduce th e question to a finite- dimensional 



matrix problem. We will also, using [BJKWOOj, show that every wavelet repre- 



sentation, if it is reducible, decomposes into a finite orthogonal sum of irreducible 



representations, i.e., if the Sj operators from (3.8) are given, then there is a finite 
orthogonal splitting 

(4.1) e\z) = Y!^np 

p 

such that each of the subspaces Tip reduces the representation, each of the restricted 
representations of Ojv is irreducible, and moreover that the irreducible subrepresen- 
tations which do occur are mutually inequivalent (and therefore disjoint). It is this 
last property of inequivalence of the irreducible subrepresentation which amounts 
to the fact that the commutant of the original representation from ( [3.8[ ) is abelian. 
Let H := £'^{Z), let B{Ti.) denote the algebra of all bounded operators in H, and let 
an arbitrary wavelet representation. Then the commutant is 

(4.2) 0'n = {X e Bin) ; = XT, Vi} = |x e B{n) ; T^XT* ^ x| . 
Of course, there are many representations oi On such that the corresponding 



commutant is not abelian, see for example |DKS99|, but the abelian prop- 
erty (i.e., that the decomposition into irreducibles is multiplicity-free) , is specific to 
the wavelet representations; see Sections ^ and ^ below. The proof of the abelian 
property is based on a lemma regarding a certain matrix which turns out to be 
diagonal with respect to a basis which is a finite subset of the Fourier basis 

(4.3) {z" ; n e Z} (also denoted e„ (z) := z"), 

or equivalently the canonical basis vectors in ^^(Z). This lemma in turn depends on 
a sublemma about a finite set of projections Pi, . . . ,Pg in Hilbert space Ti. Recall 
P G BCH) is a projection iff P = P* = P^. However, there are more details to the 
full argument, and they will be taken up in Sections ^ and ^ below. 

Lemma 4.1. Let Pi, . . . , Pg be projections. Suppose the operator 

(4.4) R = PgPg^l---P2PlP2---Pg-lPg 

is nonzero. Then R is a projection if and only if the Pi 's are mutually commuting. 
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Proof. It is clear that the operator R in (4.4) is a projection if the family Pi, . . . ,Pg 
consists of mutually commuting projections. We now prove the converse by induc- 
tion starting with two given projections Pi, P2 such that R := P2P1P2 is given to 
be a projection. Then the commutator S := P1P2 — P2P1 satisfies S* — —S. Using 
that R^ = R we conclude that = 0, and therefore 5 = 0; in other words, the 
two projections Pi,P2 commute. 



Suppose the lemma holds for fewer than g projections. If R is given as in (4.4) 
then 



(4.5) 

where 

(4.6) 



R 



PgTPg 



T — Pq-l ■ ■ ■ P2P1P2 ■ ■ ■ Pq-l- 



Writing the operator T in matrix form relative to the two projections Pg and 
/ - Pg, we get 



a 



(4.7) 

with To.o 
(4.8) 
and 
(4.9) 

imply To,! (^0,1) 
reduces to 



T 



R 

P-^TP 
^ ^a 



P TP^ 
^a 

P^TP^ 
a a 



(-^y )i.i=o 



R, etc. But then (|4.8|)-([4.9[) yield the conclusion: 
(P'^)o,o = (?o,o)^ + Po.iPl.o = P + Po,i(Po,i)* 

(P^)o,o < Po,o — R 



0, and therefore Tq^ = 0. As a result, the block matrix in (4.7) 



T = 



R 






9 9 



A further calculation shows that T must then itself be a projection. From the 
definition of T in (4^), and the induction hypothesis, we then conclude that the 
family {Pi}f^i is indeed commutative. □ 



Remark 4.2. In the special case when all the projections {Pijf^i are one-dimen- 
sional, i.e.. Pi = \vi) {vi\ in the Dirac notation, and \\vi\\ = 1, there is a simpler 
proof based on the Schwarz inequality, as follows: Let P in (4.4) be given to be 
a projection, i.e., R"^ = R ^ 0. We also have P = |Ai_2A2^3 • • • Ag_i^g| Pg with 
Xij :— {vi I Vj). We then conclude that |Ai_2A2,3 • • • Ag_i,g| = 1, and therefore by 
Schwarz, there are constants Q S C, |Ci| = 1, such that V2 — Ci'^i, V3 — C2V2, ■ • ■ , 
and the commutativity of the family {Pi}i=i is immediate. But in this case we 
find, in addition, that the projections all coincide. 

5. Minimality and representations 

The representations of the C* -algebra O at are used in other parts of mathematics, 
in addition to wavelet analysis. While it is known that in general the irreducible 
representations of On cannot be given a measurable labeling, see, e.g., [BrJo97a|, 
|BrJo97b| , |Cun77| , and jBJKWOOl , there are various families of OAr-represcntations 
which do admit labeling of their irreducibles, and their decomposition into sums 
of irreducibles. We show that the decomposition into sums of irreducibles occurs 



only for the special (permutative) representations |BrJo99a| which generalize those 
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derived from the Haar wavelets. When decompositions do occur, the irreducibles 
have multipUcity at most one; see Section ^ below. The basis for our analysis is 
the presence of certain finite-dimensional subspaces K, which are invariant under 
the operators S* when the representation is defined from the SiS with relations 
(^^)-(3.10). For related OAr-representations which arise in statistical mechanics, 
see |FNW92|, [FNW94 , and Section ^ below. These finite-dimensional subspaces 
have the the significance of labeling the correlations of the sites in the quantum 
spin chain model. If it is an infinite spin model on a one-dimensional lattice, then 
K, describes the correlations of spin observables ao,<Ti, . . . with those on the other 
side, . . . , (T_2, (T_i. 

We say that a representation of On in a Hilbert space 7i is a wavelet represen- 
tation ii Ti, = (T) (= l'^ (Z)) and if the corresponding operators Si are given 
by (3^) for some QMF functions {wij^g^. By ( 3.12 )-( 3T^ ) that is equivalent to 
using polynomial functions A: T ^ \Jn (C) for labeling the representations. We 
will let 'P(T,UAr(C)) be the semigroup of such polynomial loops, loops because 
they may be viewed as loops in the unitary group \Jn (C), see [PrSe86|. We will 
use the notation A (z) — {Ai,j {z))fjl}^ for the loop-group element A: T — > U {N). 
Since the Fourier expansion is finite, there is a g such that A (z) has the form 



9-1 



(5.1) 



^(z) = ^z'=AW (A^f-i^^^O) 



k=0 



where yl'*^' G B (C^) for k — 0, . . . , g — 1. The factorization in |BrJoOO, Lemma 
3.3] motivates the name genus for g. 



Lemma 5.1. If A{z) is a general polynomial of z with values in S (C^) of the 
form (5.1), the following four conditions (5.2)-(pTq) are equivalent: 



(5.2) A (z) A (z) = Iat, z G T, i.e., A takes values in U (N); 



(5.3) ^(fe) 



with the convention that 



ifn^ 0, 

^0 i/neZ\{0}, 

= i/m ^ {0, 1,.. .,.g- 1}; 

(5.4) there are projections Pi, . . . , in B (C^), positive integers ri, . . . , 
and a unitary W G V (N) such that 

A {z) = (n;=i (Iw - + ^''^^. )) W; 



(5.5) 



there are projections Qo, Qi, . . . , Qg-2 and a unitary y G U (N) such 
that 
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S-2 

j=o 

=F^(l^^-Qo)••• 
J=o 

• ■ ■ (iw — Qj-i) Qj (iw ^ Qj+i] 

■ • • (Iat — Qg-2) , 



S-2 



Proof. We refer the reader to [BrJoOO, Proposition 3.2] 



□ 



Remark 5. 2. The case g = 2 ~ N includes the fam ily of wavelets introduced by 
Daubechies [ Dau92 | and studied further in BEJOO . Note that g — 2 yields the 
representation 

(5.6) A^°'> ^ V {In - Q) , A^^^^VQ, 



by ( |5.5| ). But then (5^) takes the form 

(5.7) = Hat-q, ^(i)*A(i) = g, 

which will be used in the Sections ^ and || below. 

In the general case, we will need the operators (alias matrices) R (fc, I) 
^(0*^(fe)^ and the representation ( |5.5D then yields 

(5.8) 



Rig-l,g-l) = Qg-2 ■ ■ ■ QiQoQi • • • Q, 



g-2, 



which were introduced in Lemma 4.1 above 



A loop A e 7'(T, UAr(C)) is viewed as an entire analytic matrix function, C — > 
AfAr(C), and we consider ( ^.l[ ) also as a representation for this extended (entire) 
function. The (unique) entire extension will be denoted A{z) as well. The estimates 
in the next corollary translate into a stability property for the corresponding wavelet 
filters, the significance of which will be established in Section || below. 

Corollary 5.3. // g is the genus, then we have the following estimate relative to 
the order on the positive operators on ; 

(min(l,|z|2))'' ' • Iat < A(z)M(z) < (max(l,|z|2))'' ' • 1^, 
valid for all z e C, where Ti-N is the identity matrix. 

Proof. The corollary is applied in Section || below, so we postpone its proof to 
Section ||. The argument is in Observation ^.9| , and it is based on the ordered 
factorizations (5.4)-(5.5) in Lemma pTl| . □ 
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It follows from the corollary and (3.12) that the system mo, mi, . . . , mAr_i of 
polynomials that makes up the multiresolution filter cannot have any other common 
zeroes than z = 0, i.e., if some zq & C satisfies rrii (zq) = for all «, then zq = 0. 

We now turn to some representation theory for the C*-algebra On which will be 
needed in the following sections. Some background references for this are [BrJo99a|, 
Eva80[ , | Pop92|, an d |ReWc98|. Our references for wavelets and filters are |Hor95 |, 



Pol90f , and [ |Vai93[ . 



Let P {g B (Ti)) be a projection. We say that it is co-invariant for some (fixed) 
representation {T^} .^^^ oi On if 

(5.9) 7;*P = PT*P for all i. 

Let 7i„ be the closed span of {z^" ; n = 0, 1, . . . }, and let P_ be the projection 
onto 7i_. Then (5.9) is satisfied for P_ and all wavelet representations T^^\ as 
follows from (3.8), ( 3.12| ), and the formula 

JV-l 

(5.10) tI^^* = J2^^)S*, 

where S* are the adjoints of the respective operators Sj in ( ^.14 ). Specifically, 

(5.11) {s;f) (-) = ^ E 



Lemma 5.4. Let E and P be co-invariant projections for a fixed representation 
T^^"^ . Suppose E < P < P-, and further that for some r e N, 

-k 



(5.12) 

Finally assume that 
(5.13) 

Then we have the following identities 



PH = span {z-^ ; < fc < r} 



ET^'^^*P 



for alii = 0,1,..., N - 1. 



N-l 



(5.14) 



S*ES,P 



jEAi,kPk 



for all j, k = 0,1, N — I, where the functions Ai_j are the matrix entries of the 
given loop, a function is identified with the corresponding multiplication operator in 
7i — (T), and Pk S^PS/^ are projections. 



Proof. It is given that both E and P satisfy (5^) relative to T^-^\ and further that 
t(^)*E = ET^^^^*P. Equivalently, by (|l|, Y^ihi^tE = Ej:,A^^,S*,P. Using 
J2i = ^j,h we get 



(5.15) 



Now multiplying through from the right with SkP on both sides in ( 5.15 ), the 
conclusion of the lemma follows. To see this, notice first from ( 5.11 )-( 5.12 ) that 

(5.16) S*kPSiP = iik^l. 



The proof of (5.16) is based on the observation that the representation S (= T^'")) 
in ( pT^ ) is pcrmutative, see lBrJo99a| . Specifically, (z^""^) = Sk,jZ-" if < 
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j < N, and n e Z, and Si (z"") 



J-nN 



z' . The desired formula (5.14) now follows 
Pk, since P is relatively co- invariant for the 
representation S = T^^") by assumption. □ 



from this and S*^PS^^P 



Remark and Terminology 5.5. The proof shows more generally that the impli- 
cation ( ^.15 ) (5.14) holds for any operator E e B (PTi) when P is specified 
as in the statement of the lemma. In B{P'H), we may then introduce the basis 



e^k.-i 
(5.17) 



and coordinates 



k.l 



If the loop A (z) is given by (^.l|), then the operators 
(5.18) R{k,l) := A''^^*A^''^ 



of Lemma 5.1 go into the calculation of the right-hand side in (5.14) as follows: 
The (r, s)-matrix entry of the matrix {S* ES^. —) X^iLo^ ^ijEA^ ^P is given by the 
following matrix product: 



(5.19) 



Xp,qR {p-r,q- s)^^ 



P,<1 

p>r, q>s 



again with the convention that the summation indices restrict to the range where 
the terms in the sum are defined and nonvanishing. 

Definition 5.6. Let {Ti}^J^^ be a representation of On on a Hilbert space Tl, and 
let /C be a finite-dimensional subspace which satisfies 

(5.20) T*IC C /C for all i. 
Hence the projection P onto /C satisfies 

(5.21) PTi = PT,P for all i. 

We say that IC is cyclic if it is cyclic for the OAr-representation, i.e., if 

(5.22) Y r,,r,, •••r,„/c = H. 



For the wavelet representations, Ti. = L? (T), the Fourier basis {2 



; n 



E Z} has 



the following property: There is an ro G N such that, for all n e Z, there is a p e N 
satisfying 

(5.23) • • • T*^T*^ {z") e span ; < fc < tq} 

for all multi- indices ii, . . . ,iq and q > p. We showed [ Br JoOO| that tq may be taken 

9N-I 



(5.24) 



ro 



TV- 1 



where g is the genus, N is t he scale, and [xj is the largest integer < x. We also 
showed that, whenever ( 5.23 ) holds, then 

(5.25) /C := spanjz^'' ; < fc < ro} 

is cyclic. It is known in general that, if some JC is minimal with respect to the two 
properties, (5.2C) and OAr-cyclicity, then 



(5.26) 



B{IC) 



|x e s (JC) ; *p = x} 
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is an algebra. The set ( 5.26| ) is the fixed-point set for the completely positive map 
(5.27) '^ltH-) = J2^A-)V*, where V, PT^^^ . 



We further showed in [ |BJKW00t that T(^) is irreducible if and only if B (JC)" 
C Ik:- In general, this set is not an algebra, but the above minimality on /C forces 



it to be an algebra, see [DKS99 



We shall need, in the later proofs, the following two results from [BJKWOOj. We 
include the statements here since they seem not to be well known in the wavelet 
community. Let tt be a representation of Ojv on a Hilbert space Ti.. 



Theorem 5.7. |BJKWOC, Section 6] There is a positive norm-preserving linear 
isomorphism between the commutant algebra 

(5.28) TT (OnY = {AeB(n); At: {x) = n (x) A for all x £ On} 
and the fixed-point set 

(5.29) B {ICf ^ {A e B {K.) ; a (A) ^ A} 
given by 

(5.30) TT (On)' 3 a — > PAP, 

where P is the projection ofTi onto K,. In particular, tt is irreducible if and only if 
a is ergodic (where a is the mapping B (JC) —^ B (JC) defined in (5.27)). 

More generally, if ICi, K.2 (with corresponding projections P^^^ and P*^^-*) are 
T* -invariant cyclic subspaces for two representations tti, 7I"2 of On on Hi, H2, and 

(5.31) =pO)7r^.(s,)k, 
for j ^ 1,2, i = 0, . . . , N — I, define p on B (/Ci, /C2) by 
(5.32) 



p{A) = Tv}'^AVy^' 



i 

Then there is an isometric linear isomorphism between the set of intertwiners 

(5.33) {A e 6(^1,^2) ; ATTiix) ^t:2{x) A for all x e On} 
and the fixed-point set 

(5.34) {BeB{ICi,IC2);p{B)^B} 
given by 

(5.35) A — >P = P(2)AP(1). 



Theorem 5.8. pJKWOOj , Theorem 3.5] Let (fi ^ ■ V* be a normal unital 

completely positive map of B (/C) . Then 

{¥„¥*}' ciB{nr . 

Furthermore, the space B (Ti.)''' contains a largest *-subalgebra, and this algebra is 
{Vr,V*}'. 

We are now ready to state and prove the main result of the present section. Its 
significance becomes more clear when it is seen in the light of the two previous 
Theorems 5.7 and In particular, we will show in Section below that Theorem 



5.8 is applicable in verifying irreducibility, as we will show that B (JC) is generically 
an algebra for the wavelet representations. 
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Theorem 5.9. Let T^^^ he a wavelet representation of On on TC — (T), and 
assume the genus of A is g. Let tq be as in (5.24), and let P be the projection onto 
JC :— span {^z~^ ; < fc < tq} . Suppose there is a second projection E € B (JC) such 

that ^ E ^ P, and E commutes with T^'^^ * P for all i = 0,1, N — 1. Then 
it follows that E is diagonal with respect to the basis {z^^ ; fc = 0, 1, . . . , ro} in K.. 
Moreover, A (z) has a matrix corner of the form 



(5.36) 







V 







\ • • • / 

where V G JJm (C), and where the exponents Ui of the diagonal corner in A{z) 
satisfy < < 5 — 1 for all i — 0,1, M — I. 

Remark 5.10. The loops A: T — > Uat (C) which do admit nontrivial projections E 
as in the statement of the theorem are described in detail in Definition 6.6 below, 
to which we refer. Hence, the existence of such a projection E means that it is 
possible to "split off" a matrix block in A{z) which is in diagonal form. 



5.9. Let V, := PT} 



(A). 



Proof of Theorem 
V*E, or equivalently ET^ 
we have 

(5.37) {S*ES,)^^^ = ^ 

p^r, q>s 

The J, fc-indices are in {0,1, 
identities 



P ^P 



(A). 



Suppose E e B (PH) satisfies EV* 
E. Then by Lemma ^.4| and Remark 5.5 



P,qR {p-r,q- S)fe^j- = Xp+r,q+sR (p, q)^^^ 

p.q 
p>0, g>0 



, TV — 1}. For the matrices R{p,q), we have the 



(5.38) 



Yr{p,p)^1n and ^i?(p,p + = if / 7^ 0. 

p 



See Lemma |5.l| above. The terms on the left-hand side in (5.37) are 



(5.39) 



{SjESf, 



again with the convention that the terms are defined to be zero when the subscript 
indices are not in the prescribed range. 

If 7^ 0, consider the lexicographic order on the subscript indices of the cor- 
responding matrix entries Xp^q (in (5.17)). The range on both indices p,q is 
{0, 1, 2, . . . , To} where tq is determined as in Lemma 5.4, see also (5.24). Then 
pick the last (relative to lexicographic order) nonzero Xr^s, i.e., r, s are determined 
such that 



(5.40) 



X 



p+r,q+s 







if p > or g > 0. 



It follows that there are only the following possibilities for this (r, s): 



(0,0) 

(1,1) 
(2,2) 



E 
E 
E : 



X, 



0,0 \z") {z' 



01 



Xi,i|z-i)(z-i|, 

Xl,l|z-l)(z-l|+X2,2h-')(^' 
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If (r, s) = (0,0), then, using ( f3.37 ) and ( |3.39 ), we arrive at the matrix identity 



X, 



0,0 



(5.41) Xo.oi?(0,0) 



V 







e Mn (C) 



where X^ q ^ 0, 



/ 



and therefore _R(0, 0) = |eo) (£o| where Eq is the first canonical basis vector in C^. 
By Lemma 4.1, (5.8), and Remark 4.2, we conclude that Qj- > \eo) (eqI for all i, 
and therefore 



/ 1 



(5.42) 



A{z)^V 



V 







Biz) 



J 



for some V e Vn (C) and B eV {T, \Jn-i (C)); see Lemma 5.1. 

If (r,s) — (1,1), then, using again (5.37) and ( 5.39| ), we arrive at the matrix 
identity 



/ 









eMAr(C), where Xi4=^0, 







(5.43) Xi4i?(0,0): 

and therefore i? (0, 0) 
4.2, we conclude that > |eAr_i) (e^v-il for all i, and therefore 



V 



\ 





^1,1 / 



^N-i) i^N-il- Using again Lemma 4.1, (5.8), and Remark 



(5.44) 



A{z)^V 



C{z) 



V 











1 / 



for some y e Uw (C) and C G 7^ (T, U^v-i (C)). 

The reason for why the matrix of E has diagonal form relative to the natural 
Fourier basis is as follows: Let N = 2, for simplicity. (The argument is the same, 
mutatis mutandis, in the general case.) Then pick the last term (r, s), r ^ s, with 
Xr,s 7^ 0, where again "last" refers to the lexicographic order of the matrix-entry 
indices, see (5.40). We then get, using (5.37) and (5.39), the following matrix- 
identity (where we specialize to (r, s) = (0, 1)): 



^o,ii?(0,0) 




^0,1 



e A'h (C) 



and Xq^i ^ 0, as mentioned. This forces i?(0,0) = (lo)) which is impossible by 
Lemma 5.1, since i?(0,0) is positive and (? q) 

Let N > 2, and suppose (r, s) = (2,2), i.e., assume that X2.2 7^ 0, and 
X2+p,2+q = if p > or g > 0, referring to the lexicographic order. Then by 
the same argument which we used in the earlier cases, 



(5.45) 



X2,2(i?(0,0)),,^. = iX2N^,,2N-k)^,^^o 
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But all the double indices {2N — j, 2N — k) of the matrix on the right are strictly 
bigger than (2, 2) in the lexicographic order, and we conclude that 



(5.46) 



R (0, 0) = in A'l 



.'IN 



The formula for R (0, 0) then yields QqQ^^ ■ ■ ■ Q^-i = 0. Moreover, the additional 
restrictions are: 

Xo,oi? (0, 0) + XisR (1, 1) + X2,2R (2, 2) = Xo,o ko) (£o| € Mm (C) 



as in (5.41), and 





( " ■ 


• 








Xi,ii?(0,0)+X2,2i?(l,l) = 


• 


• 










• 


• 


X2,2 







^0 . 


• 





^1,1 



Now substituting R (0, 0) — 0, we arrive at 

(5.47) (1,1)+ X2,2R (2, 2) = Xo,o |£o) {eo\ 
and 

(5.48) X2.2R (1, 1) - X2,2 \sN-2) {eN-2\ + \eN-i) (ew-i| . 

Since E = Xq.o |1) (1| + Xi^i \z^^) {^^'^\ + ^2,2 (-^^^l is a projection, and 

X2.2 7^ 0, we must have ^2,2 ~ 1 and Xq.o and Xi i € {0 , Ij. The conclusion of 
the theorem can then be checked case by case, using ( |5^ ) and (|1|). 

In general, let v'C^^s be the last (in lexicographic order) nonzero term, and 
assume s > 2. Then by ( 5.3^ )-( 5.39 ), we get Xs^s-R(0, 0) = 0, and therefore 
R (0, 0) = as before. Using this, the equation for the (s — 1, s — 1) term is then 

-^s,s^(l,l) — {X{^s-l)N-j,{s-l)N-k) j k^Q ■ 

If (s — 1) {N — 1) < N , then all the entries in the matrix on the right must vanish, 
and we get R (1, 1) = 0. If not, we proceed as in ( 5.48 ). If R (1, 1) ~ 0, we go to 
the (s — 2, s — 2) term, viz.. 



(5.49) 



X,,,i?(2,2) = (X(,_2) 



2)N-j,(s^2)N-k)jj^^Q ■ 

Eventually the matrix on the right will have nonzero terms, starting with Xs^s, and 
terms before that in the lexicographic order. Suppose, for example, that the matrix 
on the right in (5.49) has nonzero entries. Then the equation for the (s — 3, s — 3) 
term is 



Xs-us-iR (2, 2) + Xs^sR (3, 3) = {X, 



{s~3)N~j,{s~-3)N~k) , 

and the argument is done by a case-by-case check, using that the coordinates 
Xo,o,-'^i,i, . . • are in {0, 1} while Xg^s = 1. 

There is a similar argument, based on the reversed lexicographic order, starting 
with (N — I, N — 1), which will account for a possible lower right matrix corner of 
diagonal form. This completes the proof of the theorem. □ 
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Remark 5.11. (Permutative Representations) The form 



(5.50) 



A{z) = V 







\ • • • / 

V e \Jn (C), in the conclusion of Theorem |5.9| corresponds to the representations 
of On which permute the basis ele ments {z" ; n E Z} ioi H ^ (T); they are 
studied more generally in [BrJo99a| under the name permutative representations. 
We also met them, in a special case, in Remark 2.7 above, in connection with the 
"stretched out" Haar wavelets. So the conclusion of Theorem 5.£ is that the wavelet 
representations which are not of this form are irreducible. 

Now for the details: Let T^"^^ be the representation of On corresponding to A (z) 
in ( ^.5C| ). The element V eVn (C) defines an automorphism of Oat, denoted ay or 
Ad {V). Let D (z) = V-^A (z) be the diagonal factor in ( |53o| ). If tt^^) (s,) = t/^^ 
and 7r(^) (s,) = xj;^^ are the corresponding representations, then it follows that 



(5.51) 

which means that tt^^^ and tt*^^^ have the same decomposition into sums of irre- 
ducibles, corresponding to irreducible subspaces of (T). The formulas for the 
operators Tj^'^^ are as follows: 



(5.52) {z') 



— ^ 1 



fc e z, i = 0, . 



,iv~i. 



which justifies the "permutative" label; in other words, both the operators 
and their adjoints permute the basis elements of the Fourier basis {z*^ ; fc G Z} 
for (T) . The de composit ion structure of these representations is worked out in 
BrJo994 ; see also ^KS9|. 



Remark 5.12. Note that if A'' > 2, the n som e representation T^^^ may be re- 
ducible even if A is not itself of the form ( 5.50| ); it may only have a matrix corner 
of this form. Take, for example. 







1 1 








(5.53) 


A{z)^ 




I 


1/V2 
z/V2 


z/V2 



i.e., = 3, 5 = 3. Then T^-^^ is a reducible representation of O3 acting on (T), 
and, in fact, the Hardy subspace H'^ C (T) reduces this representation. For 
more details, see Section |6| below. 



(6.1) 



6. Irreducibility 

In this section we consider as in Lemma ^.4| ( ^.25| ) the finite-dimensional subspace 

.gA^ - 1 



K := span{z"'= ; < fc < ro} C i-'(T), 



N - 1 



defined from a polynomial loop A(z) of scale size A^ and genus g, and we show 
that the irreducibility property for the corresponding representation T'^'^^ of Oat is 
generic, i.e., it holds for all A except for a subvariety of smaller dimension, once 
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N and g are fixed. In order to apply the results in Sections ^ and some more 
details are needed regarding the subspace /C, and they are taken up in Section |8[ 
We begin with some notation and a lemma: 

Notation 6.1. Let e„(z) :— z^\n e Z, denote the Fourier basis in L'^{T) 
finite subsets J C Z, set (J) := spanjej ; j G J} C L^(T). If 



For 



(6.2) 



Jo = {0,-1,-2, 



set /c := (Jo). If r = r(^) is a wavelet representation and ro is as above, we note 
fBrJoOO , Proposition 5.5] that /C is cyclic. 



Lemma 6.2. Let A G V{T,\J n{C)) be a (polynomial) loop of genus g. Then 



(6.3) /C=({0,-1, 



-ro 



}) = span {z '^;0<A:<ro}, ro = 



gN-i 

N - 1 



contains no one- dimensional subspace which is both T^^^ * -invariant, and also cyclic 
for the representation of On on L^(T). 

Proof. To show that a subspace /C is minimal in the sense specified in the lemma, 
we must check that whenever 



ICo ^ K. 



r/^^*/Co c/Co 



for i = 0, 1, . . .,7V- 1, 



(6.4) 

is a subspace satisfying 
(6.5) 

then /Co cannot be cyclic for the representation T^"^^ of Oat, i.e., it generates a 
cyclic subspace which is a proper subspace of L^(T). The cyclic subspace generated 
by /Co is the closed subspace spanned by 

(6.6) T^[^^ ■■■T^^flCo forn = 0,l,..., 

' ,i' 



and all multi-indices (ii, . . . ,i„). This follows from ( |6.5| ), and we will denote this 
space [Oat/Co] . We will prove the assertion by checking that if ( |6.4| )-(3.5) hold, 
then there is an mo G Z such that [Oat /Co] is contained in the closed span of 
{z'' ; k > Too}. Note that this integer too might be negative, and also that [©Af/Co] 
might well be a proper subspace. 



Now let /Co be given subject to conditions (6.4)-(6.5), and suppose (as in the 

= 1. Then by (6.5) there are Ai S C with 



lemma) that dim /Co = 1. Let ^ G /Co 
(6.7) 
or equivalently, 

(6.8) e(^) = VA-TO|^)(z)e(z^) 



Using |Jor99a , |BrJo97b| we conclude that 

(6.9) ^(z) = z^'' for some k, 

after adjusting with a constant multiple, and moreover that 



(6.10) 



(A) 



(^) 



.fc(JV-i) 
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Setting a{z) := (l, z, . . 
(6.11) 

Now pick j e {0, 1, . . . 



, this may be rewritten as 
— 1} such that — fc = j mod A^, and apply the operators 



such that 
(6.12) 



iV — 1, to both sides in (S.ll). It foUows that there is some rii e 



(A I A{z)ej) = and (A | A{z)ei) =0 iil^j. 



Since |(A|^(2;)£j)| < ||A|| || A(z)£j | = 1, the first part of (6.12) imphes equahty in 



a Schwarz inequahty. Then ( p.l 
(6.13) 



yields A{z)*X 
A{z)ej ^ z™A. 



or equivalently, 



Using the formula in Lemma 5.1 for the coefficients in A{z), and Lemma 4.1, we 
note that ( 6.13| ) implies 



(6.14) 



3' 



and in particular A — Vsj, where V G \Jn{C) is as in Lemma 5.1. 

For the convenience of the reader, we sketch the argument, but only in the 
simplest case m = 0. With the notation of Lemma we get 

(6.15) A^°hj = A, A^'^hj =0, l<k<g, 

where j is the (fixed) number determined from ( |6.11 ) as described. Introducing 
the projections Qq, Qi, • • • e B{C'^) from Lemma 5.1, the first part of ( |6.15 ) then 
reads 



(6.16) 
and since 

(6.17) 



VQiQi---Qj_,e,=X, 



\M\ = 




= 1, 



IQo *3i ' ' ' Qn-i^j II — 1- 



■ Qg^iBj is in the range of Qq, and Qq{Qi 



Hence Qf- 
We get 

(6.18) ||Q^...Q^_^£^.|| = 
and by induction, 

(6.19) Qf-i£j = Q^-2^j = . . . = Q^Ej - >^o - -^j 

and therefore Qk£j = 0, and A{z)ej = Vej = A. This proves the claim ( 6.14 ). 
Since we wish to prove that [OnICq] is contained in z™Ti.^ for some m € Z, where 
H+ : = span {z^ ; < fc} is the Hardy space in L^(T), we may assume that m in 
( 3.14 ) is taken as m = 0. Since the invariant subspaces for a representation tt of Ojv 
are the same as for tt o av where ay — Ad V, we may replace A{z ) with V^^A{z), 
or equivalently, reduce to the special case A{z)ej = ej of formula ( |6.14 ). Since the 
matrix of the basis permutation Eq ^ ej is in UAr(C), the same argument leaves us 
with the simpler case A{z)eo = Eqi or equivalently. 



/ 1 



(6.20) 



A{z) = 



V 











AN-isiz) 



Ai,Ar_i(z) 
An-1.N-i{z) j 
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Then [ BrJoOO| , Theorem 6.2] imphes £,{z) — z° = 1, or equivalently k = in (6.9), 
and so [Oat^] C H+. Since /Cq = C^, we have proved that /Co is not cychc, i.e., the 
representation on the single vector ^ does not generate all of (T) . This concludes 
the proof of the lemma. □ 



Remark 6.3. If A e V{T,\Jn{C)), then the wavelet representation T^-^) may or 
may not be irreducible, as a representation of On on L'^{T). It is not irreducible, 
for example, for the Haar wavelet. Nonetheless, we will show that when N and 
g are given, then irreducibility holds generically for T^^) as A ranges over all 
^^(T, UAr(C)), i.e., has scale number N and genus g. 

Now for r('^) irreducible, of course every ^ e L'^{T), ^ 7^ 0, will be cyclic, so 
cyclicity will then not be an issue, but rather the question of when 1^ satisfies 

(6.21) i;^^^*^ e ce 



for all i. 



The lemma states that ( |6.21 ) cannot hold for ^ 7^ if T^^^ is irreducible. The issue 
is then instead to find minimal subspaccs /Crcd such that 



.22) 



rp{A)* ^^^^^^ (- /Crod- 



These subspaces are relevant for the algorithms which are used in the construction 
of and the analysis of wavelets, as we sketched in Section |l|, as the subspaces /Crod 
in (T) = £2 correspond to subspaces in the associated multiresolution subspaces 
in L2(M) (See Table |l|). We first addressed ( |6.22|) in | |BrJoOC| |, but the minimality 
was not considered there. We also note that ( |6.2^ ) has applications for different 
representations of On, and is there connected with finitely correlated states in 
statistical mechanics, see prJo97a[ |FNW94|, p^"NW92| . The minimality issue was 



also considered in [ pKS9!i | in a different context. We noted in ( 5.25| ) that /C 



span{2;°, z"^, . . . ,z~^°}, with ro as in (5.24), satisfies (S.22). Let P or Pjc denote 
the projection onto the subspace /C. We will consider subspaces of /C which still 
satisfy ( |6.22 ), are cyclic, and minimal with respect to (6.22) and cyclicity. If, for 
example, g = 3 and = 2, we will show that for some A e 7'3(T, U2(C)) we may 
have 



(6.23) 



/Crcd = span 



This is a little surprising since then tq = 5, and so JC in (5.25) is 6-dimensional 



Corollary 6.4. Let Jo be as above in ( |6.2| ), and consider the two finite- dimensional 
subspaces ICq — (^o\{0}) and fCi — {Ja\{^'f{)\) in Lp'iT). Then JCq is non-cyclic if 

(6.24) Tl;^^*eo G Ceo for all i. 

Suppose N — 1 divides gN — 1. Then JCi is non-cyclic if 



(6.25) 



T, 



(e-ro) e 



for all i. 



Moreover, ICq is cyclic if Xq = J2i 



AO) 



= 0. 



Proof. The two vectors cq and e_ro, corresponding to the endpoints in Jo, are 
special in that 



(6.26) 



PT^^^^eo e Ceo, 
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and when iV — 1 divides gN — 1, 



(6.27) 



e Ce_ 



for all i. 



It follows that, when (6.24) holds, then 

(6.28) [On ({0})] ® [O^ (Jo\{0})] - L'(T) 

where ({0}) = Cz" = the one-dimensional space of the cons tants , and neither 
of the two subspaces in this orthogonal sum is zero. Hence ( |6.24 ) implies that 
^0 — ('A)\{0}) is not cyclic. The same argument proves that /Ci = (•A)\{~^o}) is 
not cyclic if ( |6.25 ) holds. 

We also note that ( |6.24| ) holds if and only if 



(6.29) 

or equivalently, 
(6.30) 



= 1, 



i.e., Ao [A) = 1, 



1. 



(The issue is resumed in Remark 6.5 below.) 

We now turn to the converse implications in the corollary, doing the details only 
for /Co = (Jo\{0}) . If (|6^) does not hold, then 



(6.31) 



(0) 



satisfies Ao < 1. Using the following formula. 



(6.32) 



eo 



T/^)(e_,)) = (T, 



eo 



(z'=iAo(z))=4:;^ 



z-'=) = ^A,,o(z) 
(fe) 



we therefore have the conclusion: For each z, T^^^* is in /C, and it splits according 
to the sum K, — Ceo + /Co a-s follows: 



i;(^)*eo = AfJeo + 6, 0<i<iV, 



(6.33) 

where S^i £ /Cp is computed according to formula ( |6.32 ). Applying PT^^^ to both 
sides in ( 6.33 ), we conclude that 



eo 



or equivalently, 

(6.34) (1 - Ao) eo e P [Oat/Co] . 

Since in the second part of the corollary, Aq = by assumption, we conclude from 
( |6.33| ) that eo G [Oat /Co], and the inclusion 

(6.35) ({0,-l,...,-ro}) c [Oat/Co] 



-Jo 



follows. Finally, we get 



(6.36) 



L\T) = [On (Jo)] C [Oat/Co] C L\T), 
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which proves that the reduced space /Co is then cychc, and the coroUary is estab- 
hshed. □ 

Remark 6.5. To summarize, T^"*^ is given by some A e V{T,\J n {€■)), where the 
coefficients A'^°\A^'^\ • • • e Mn{C) satisfy A^^-i) ^ q and the conditions in (O) 



and Lemma 5.1. The conchisions about cychcity in the previous corollary may be 



restated as follows in terms of these matrices: 

(6.37) eo ^ [On ( Jo\{0})] <=> A has the form 



* 1 




1 




1 






1 




* 1 




i 




i 






i 






M 




/ 






c 




L 


* 1 




1 




1 






1 





















Ad 



A(9- 



i.e., all in the first of the first columns; and 



(6.38) 



[On ( Jo\{-ro})] 



A has the form 





1 




1 




1 






1 




1 * 




1 




1 




1 






1 




1 * 


M 


1 


/ 


1 


s 


1 




c 


1 


L 


1 !); 























Ad) 



Ai9- 



i.e., all in the last of the last columns. It follows from Theorem 5.9 that the 
representation T''^^ is reducible if either one of the conditions ( 6.37 ) or ( 6.3^ ) holds. 

We can show, using [BrJoOO, Theorem 6.2], that a wavelet representation T'"^^ 
satisfies (3.37) if and only if there are V € UAr(C) and B g ■p(T, UAr_i(C)) such 
that A ha s the form (5.42). There is a similar conclusion concerning the other 
condition ( 6.38| ). 

Let us say that, for fixed N and g, a property is generic if it holds for all 
loops A{z) of scale N and genus g, except for A in a variety of lower dimen- 
sion. Then we conclude from ( 6.37 )^( 3.38 ) that /Crcd = ({"li • ■ ■ : ^(''o — 1)}) or 
({—1, . . . ,—{ro — 1), — ro}) is cy clic fo r a generic set of loops, when g and N are 
fixed. The process described in (6.31) of elimination starting with the elimination 
of, if possible, and — rp from ({0,-1, . . . , — rp}) to get a smaller space, say /Cq 
such that 



(6.39) 
and 



i;*^^*(/Co) C /Co foralH, 



(6.40) /Co is OAT-cyclic in i^(T), 

may be continued, subject to certain spectral conditions on the given loop A. These 
conditions are generic in the same sense; for example, as noted in ( |6.31 ), eo can 
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(6.41) 



be eliminated (so that the smaller /Cq will still satisfy ( 6.39D ~( |6.40 )) if and only 
if i?(0, 0)0.0 < 1- There is a similar spectral condition for the elimination of two 
vecto rs eo and e_i, i.e., for getting /Co = ({—2, —3, . . . , — ro}) to also satisfy ( |6.39 )- 
( |6l0| ): For the iV = 2 case, this condition is that the 2-by-2 matrix 

i?(l, 1)0,0 i?(0, 1)0,1 
i?(l, 0)1,0 i?(0,0)i,i 

does not have 1 in its spectrum (recall R{k,l) :— A^'^*^^'^)). The argument is the 
same as before, even if A'^ > 2, mutatis mutandis. If 1 is not in the spectrum, and if 
i?(0, 0)0,0 = 0, then we show that both the vectors eo and e_i are in the cut-down 
of the cyclic space, i.e., in 

(6.42) P[02({-2,-3,...,-ro})], 

which is then rj*-"*''*-invariant. (Here we use the symbol P for the projection onto the 



Ng-l 
N-l 



IfiV = 2, then. 



subspace /C spanned by{z '^;fc = 0,l,...,ro} where ro 

of course, JC is of dimension 2g.) Hence, the smaller subs pa ce ({ —2, —3, . 
({0,-1, —2, . . . , — ro}) will also satisfy conditions ( 6.39| )-( |6!40 ). 

To illustrate the spectral condition more explicitly, we need g > 2. In the case 
g — 3, N = 2, there are V G U2(C), and projections P, Q in C^, such that 

(6.43) A{z) = V{P^+zP){Q^ + zQ), 

and we get 

i?(o,o) = Q^p^g^, i?(o,i) = gp^Q^, i?(i,o) = Q^p^g, 

R{l,l) = QP^Q + Q^PQ^, R{2,2)^QPQ. 
Hence, R{0, 0)o,o — holds if and only if P^Q^Eq = 0, or equivalently, 
(6.45) PQeo=Peo + Qeo-eo. 



(6.44) 



The entries of the matrix (3.41) are then 
Since 



{eo\P^Q^ei) 

1 
1 



(6.47) i?(0,0) + i?(l,l) + i?(2,2) = 

the restriction i?(0,0)o,o = therefore implies that 

(6.48) i?(l, 1)0,0 = 1 - i?(2, 2)0,0 = 1 - WPQeof ■ 



Using this, we get that 1 is in the spectrum of (3.41), so reduced, if and only if 



l-\\P^Q^e,\\') = \{eo\P^Q^ei)\\ 



(6.49) \\PQeo\ 

To solve this, let for example P and Q be the respective projections onto 

cos 9 \ , / cosp 



(6.50) 



smf 



and 



smp 



and solve for 9 and p. For these examples (i.e., in the complementary region of the 
(0, /9)-plane), we will then have /Cred :=({— 2,— 3}) satisfy the covariance condition 
as well as the cyclicity. All of the cases iV = 2, 5 = 3, will be taken up again in the 
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Appendix, where the algorithmic properties of (1.7) for the scahng function Lp are 
displayed in detail. This is an iteration based on (1.7), and the regularity of the 
corresponding x i— s- (ysg.p (x) turns out to depend on the spectral properties of the 
operators in (6.44). 

We now turn to the distinction between the diagonal elements A{z) in 
■p(T, U7v(C)), and the non-diagonal ones. We say that A is diagonal if it maps 
into the diagonal matrices in UAr(C), except for a constant factor, i.e., if there is 
some V G Uat (C) , and ng , . . . , njv_ i > 



0, such that A has the form ( 5.5C| ). The vari 



ety of these diagonal loops will be called T^diagC^T, Ujv(C)). Note that this definition 
includes 



(6.51) 

m7'diag(T, U2(C)). 



Definition 6.6. We say that a loop A E 7'(T, Ujv(C)) is purely non-diagonal if 
there is not a decomposition N = do + b + di with dg > 0, or di > 0, diagonal 
elements D,{z) e 7'diag(T, Ud, (C)), i = 0,1, B{z) e V{T,VbiT)), and V G Ujv(C) 
such that 



(6.52) 



A{z) = V 






B{z) 
D^iz) 



Easy examples of loop matrices for N — 2 and g — 3 which are not diagonal, 
i.e., do not have the representation (5.50) or ( |6.52|) for any V, are 



(6.53) 



1 



z 



and 



1 

71 



z 



1 

Both have \o{A) = 1/2. Both correspond to Haar wavelets, and both are ex- 
ceptional cases in the wider family of the Appendix below. Both have irreducible 
wavelet representations, by the next theorem. For more de tails about matr ix factor - 
izations in the loop groups, we refer the reader to | PrSc86| , and the paper [ AlPe99 |. 

We now turn to our first explicit result about minimal subspaces £ C /C, i.e., 
subspaces C which are r*-invariant, cyclic, and which d o n ot contain proper T*- 
invariant subspaces which are also cyclic. In Theorem 8.2 below, we shall then 



further give a formula for the (unique) minimal such space C. We stress that these 
results are special for the wavelet representations, and that they do not hold for 
other kinds of representations of On- 

Theorem 6.7. (a) Let A G V{T,\J NiC)) be given. Suppose it is purely non- 
diagonal, and let T'"*' be the corresponding wavelet representation of On on L'^{T). 
Then it follows that T^"^-* is irreducible. 

(b) Let rg be as in (5.24). Let optimal numbers p, q, Q < p < q < be determined 
by the spectral condition in Remark |6.5| such that 

(6.54) /C,ed = ({-P,-(P+1),...,-'Z}) 
is T^^''* -invariant for all i, and further satisfies 

(6.55) ({0, -1, . . . , + 1, -9 - 1, . . . , -ro}> C [Ow/Cred] ■ 
Then the following three properties hold: 
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(i) 7;''^^*(/Crcd) C /Crcd for all i, 

(ii) /Crcd is cyclic (/orL^(T)), 

(iii) /Crcd minimal with respect to properties (i) -(ii) . 

(c) The minimal space /Crcd from (b) is reduced from the right if N ~ 1 divides 
gN — 1, where g is the genus, and if not, it is {{—p, ■ ■ ■ , — ro}); so it is only "trun- 
cated" at one end when N ^ 1 does not divide gN — 1. 

Proof. Once /Crcd has been chosen as in the statement (b) of the theorem, the three 
properties (i)-(iii) foUow from Theorem 5.8 and ^.9| . The significance of (i)-(iii) is 
that they imply that if 



(6.56) 



a(.) 



{A)* 



then the fixed-point set B{K.rcdY ^^.ct an algebra. This is a result of Davidson 
et al. [DKS9£]. Using Theorem 5.8, we conclude that the projections in B{K.^cdY 
are characterized by the condition of Lemma 5.4. Now, by | DKS99[ , there are 
projections Ej G B{K,-ccd) such that, for each i, j, we have the covariance properties 

(6.57) 



e,v:e,^v:e„ 



where V* = Tj^^"^" P/c^^^, or equivalently. 



(6.58) 

and in addition, we have 



1^; = Pr_T, 



(6.59) 



and each subspace [ONEjICrcd] irreducible, in the sense that each [Cat E'j /Crcd] 
duces the represen tatio n On t o one which is irreducible on the subspace. 
It follows from ( |6.57| )-( 6.59 ) that the complementary projection 



re- 



(6.60) 



E, 



then also satisfies (3.57), and so in particular Ej must commute with each Vi 



Then by Theorems 5.8 



and 



5.9 



we conclude that each Ej has a 



matrix which is diagonal with respect to the Fourier basis {z '^j. Since the loop 
A{z) is picked to be purely non-diagonal, we finally conclude that the decomposition 
{Ej} of (3.59) can only have one term, and the proof is concluded. □ 



Remark 6.8. Even if the assumption in Theorem 6.7, to the effect that A be 
purely non-diagonal, is removed, we have the decomposition into irreducibles, and 
these irreducibles [0]\[EjlCrcd\ are mutually disjoint, i.e., inequivalent representa- 
tions when j ^ j' for two possible terms j, j' in a decomposition. This follows 
from the Theorems ^.8| and 5.9, which state that the projections Ej are all diagonal 
relative to the same basis (see also Theorem 8^ below!). So in particular, B[lCicdY 
is abelian when /Cred is chosen subject to conditions (i)-(iii) in the statement of 
Theorem 6.7. This means that the corresponding decomposition of T*-"^-* into a sum 
of irreducible representations of On is multiplicity-free. 
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Example 6.9 (An Application). Even though we hst only the scaling functions 
Lp (x) in the examples in the Appendix, the wavelet generator ip (x) is significant. 
But it is not unique: We can have a loop A of genus 2, and a different one B 
of genus 3, which have the same (p. Then, of course, there are different wavelet 
generators, say ■0^ and ip^. To see this, take (p as follows (see also Remark 2.7): 



(6.61) 



A loop A in diagonal form giving this ip is 



V2 I 1 -1 



1 

z 



genus g = 2. 



This is of the form (5.50). The corresponding wavelet generator ■0^ is then 



(6.62) 




But setting 



(6.63) 



then this loop has the same (p. Since 



(6.64) m[^^ (z) = — (1 + z3) = m^"^' (z) , m[^^ (z) 

v2 



z^m^^^ (z) , 
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the corresponding wavelet generator jp-^ is now different from ■0^ only by a trans- 
lation. 



(6.65) 



In fact, ^pg{x) — '0^ (a; — 1). However, the most striking contrast between the 
two loops A and B is that the minimality question comes out differently from one 
to the next: The representation T''^) of O2 on (T) is reducible, while T*^^^ is 
irreducible, i.e., there are no nonzero closed subspaces of (T), other than (T), 



which are invariant under all T^'^'^ and x}^"'* . (Or, stated equivalently, by (4.2) 



we have the implication Y.^T■^'' XT!>"^* =X,X d B{L'^iJ)), ^ X e C1l2(t).) 



-n(B)* _ 



The two conclusions for T^^^ and T*^-^^ follow from Lemma |672| , Corollary 6.4, and 
Theorems 5.9 and |6.7|, respectively; but Theorem |8.2| is also used. What is perhaps 



more surprising is that the matrix loop 



(6.66) 



A{z) Bi 





/ 1 


z 





\ 


I 


z 


















1 


z 




V 





z 





(see ( |6.63D ) in U4 (C), i.e., = 4 and g 
which acts irreducibly on (T) . 



3, defines a representation r(-^) of O4 



Our general result in this paper is that the wavelet representations are irre- 
ducible, except for isolated examples of Haar type, such as tpj^ in ( 3.62 ). But ( 3.64 ) 
above shows that even for the reducible ones, irrcducibility can still be achieved, 
if only Z-translations are allowed; see ( |6.65| ). The following result is a corollary of 
Theorem 3/?, and it helps to dist inguish the wavele t representations T^'^^ from the 
more general representations of FNW9S , FNW94 associated with finitely corre- 
lated states in statistical mechanics. It is a crucial distinction, and it i s co ncerned 
with the completely positive maps a which are described in Theorems p. 7] and 5 



In | FNW94 ], the representations are determined by maps a which possess faithful 
invariant states, and these states play a role in the proofs of the results there. Our 
next corollary asserts that such faithful invariant states do not exist for the wavelet 
representations . 

Corollary 6.10. Let T'^^^ be a wavelet representation of On on (T) which sat- 



isfies the conditions in Theorem 6.7, and let cr^'' {■) — i') corre- 

sponding completely positive mapping of Theorems 5.7 and 5.S. Then there is no 



faithful state p on B (IC) which leaves cr^-* invariant, i.e., which satisfies 
(6.67) 



(-4) 

poa^ 
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Proof. We will restrict to the case N — 2, although for g ~ 2, we cover arbitrary N 
in fBrJoOOl . (I f g = 2 then JC = (eo, e_i, e_2, e-s). Setting E^k-i |e-fc) (e_i|, 
we showed in [BrJoOC] that the density matrix D given by D = AAr_2-E'-i.-i + 
(1 — Ajv_i) E-2-2 satisfies a* (D) = D, where a = ct^^ and :— R (0, 0)- ,-, and 
where a* is the adjoint oi a: B (/C) B {JC) with respect to the trace inner product. 
Defining the state p on B {JC) by 

(6.68) p (X) := trace (DX) , VXeS(/C), 

we check that p satisfies ( |6.67| ). We know from jBJKWOOll that ker (cr - 1) and 
ker {a* — 1) have the same dimension. But T'^^^ is irreducible by Theorem 3.7 



when < Ao < f . Hence ker (ct — 1) is one-dimensional by Theorem 5.7, and there 
are therefore no other states p satisfying ( |6.67 ). But the state p in ( 6.68| ) is clearly 
not faithful, and the proof is complete, in the special case g = 2.) 

We now turn to the details for N ^ 2, g ~ and it will be clear that they 
generalize to arbitrary g. li N — 2, g — i, we get JC — (eo, e_i, e_2, e_3, e_4, e-s) = 
and r/ "^^ *e _fc may easily be computed; see, e.g., the details in Section ^ below. 



especially (|8!26|)-(^). As a result, we get a* {E^u^^i) = \T*e^k) {T*e^il and 
therefore 



(6.69) 
(6.70) 
(6.71) 
(6.72) 



a* {E_2,-2) 
<y* {E-3^-3) 
a* (-B_4,-4) 



^i?(?,fc),,£;_i_fc, 

k,l 



-l-l ■ 



E 

k,l 



^i?(?,fc)ii£;_2-fc,-5 

k,l 

i? {I, fc)p g £'-2-fc,-2-i , 



-2~l ■ 



k,l 



where the fc, / summations are both over {0, 1, 2}. In addition, by ( 6.26 ) and ( 6.27 ), 

(6.73) a {Eofi) ^ XoEofi, and cr (i;_5 ^5) = Aq^^-s _5, 

where Ao = Aq {A) = i?(0, 0)Qg. So the complement of (£^0,0, ^'-5, -5) in B {JC) is 
invariant under a* , and the element D which is fixed by a* must be diagonal in 
the Fourier basis, by Theorem |]^. Using Lemma 5.1 and the argument from the 
previous step, we then check that a density matrix D may be found in the form 



.74) D = SiE^i^^i + ^2^^-2,-2 + S3E^3^^3 + 5iE^i^^4 



1, 



such that the state p{-) = trace (D • ) on B {£,) will satisfy (6.67). But if < 



Ao {A) < 1, the wavelet representation T^^'' is irreducible, and so (6.67) has no 
other state solutions. Finally, it is clear from ( |6.74 ) that p is not faithful. □ 



7. FILTRATIONS IN 'P(T,U2(C)) AS 
FACTORIZATIONS OF QUADRATURE MIRROR FILTERS 

Since V{T,\J n{C)) has multiplicative structure, it has ideals, and since the uni- 
modular polynomials, i.e., T ^ T, are monomials, we may reduce the consideration 
to the ideals z'=P(T, Ujv(C)), fc = 0, 1, 2, . . . . 

In view of the examples, we specialize the discussion to the case N — 2, but the 
arguments work generally. 
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Figure 1. B{z)eVg (T, Ua (C)) vs. m'^^^ (z) = zmf (z) € P^+i (T, U2 (C)) 



In this section, we explain how the subspace K. :~ ({0, —1, . . . , — (2g — 1)}) in 
(5.25) is reduced first to the smaller one ({—1, —2, . . . , —{2g — 2)}), and then fur- 
ther to ({— 2, . . . , — (25 — 3)}), in the case N — 2. Returning to the semigroup 
7'(T, U2(C)), we note that it has a natural filtration of ideals: 



(7.1) 



z7'(T,U2(C)) D z^r{T,V2{C)) D 



A loop A{z) is in z'=7'(T, U2(C)) if and only if there is some B{z) £ 7'(T,U2(C)) 
such that 



(7.2) 
Since 
(7.3) 



A{z) = z''B{z), zeT. 
mf\z)=Y,B.A^^)z\ 



we get m\^\z) — z'^'^ m['^\z) , and for the representations 



(7.4) 



where M^2k denotes multiplication by z^*^ on the Hilbert space (T) . Despite this 
simple relationship between T^"^^ and T^^\ the irreducibility question can come 
out differently from one to the other. 

If A{z) — zB{z), and B is of genus g, then A is of genus g + 1, but it has 
vanishing first and last columns in its representation, as is clear from Figure ^ 
Specifically, suppose ml'^^ (z) = zm[^'^ (z) for all i; then we have the following 
system of identities: 

(7.5) 

4(0) _ n - R^") 4(1) - r(°) 4(1) - r(i) 
for i = 0, 1, and so the matrix 

i?(l, 1)0.0 i?(0, 1)0,1 

i?(l, 0)1,0 -R(0, 0)1,1 



of Section ^ takes the following form: 

(7f\\ ( ^-4(1,1)0,0 i?A(0, 1)0.1 
^'■''> I i?A(l, 0)1,0 i?A(0,0)i,i 



i?B(0, 0)1,1 i?i3(0, 0)1.0 
i?B(0, 0)0,1 i?i3(0, 0)0,0 
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Moreover, a given A £ Vg+iiT, U2(C)) has the form m-'^'' (z) = zm[^^ (z) for some 
B eV{T, U2(C)) if and only if 

(7.7) Ao(A)(:=i?A(0,0)o,o) = 0. 

Putting this together, we get the foUowing result: 
Proposition 7.1. (a) Let 

AePg+i(T,U2(C)), 

and let P be the projection onto the subspace IC. Then the following three 
conditions, (i), (li), and (iii), are equivalent: 
(i) Ao(v4) = 0; 

ml [z] = zm] (zj ; 
(iii) eoG [02({-l,-2,...,-2(.g-l)})]. 
(b) The following two conditions are equivalent: 
(i) Ao(^) = 1; 



(ii) there is some 



yeU2(C), &eT, 



such that 

(7.8) A{z) = V 



1 

hz9 



for all z G T. 

(c) The following two conditions are equivalent: 
(i) Ao(^) < 1; 

(ii) 

(7.9) eoeP[02({-l,-2,...,-2g})], and 

e_(2g+i) eP[02 ({-l,-2,...,-2g})]. 

(d) Suppose Ao(^) = 0. Then the following thre e co nditions are equivalent: 



(i) 1 is not in the spectrum of the matrix (7.6); 
(ii) 

(7.10) eo,e_i e P[02 ({-2,..., -2.9 + !}>], and 

e_2s-i, e_25 e P [O2 ({-2, . . . , -2.9 + 1})] ; 



(iii) the loop B in nri'^^^ (2) ~ zm^^^ (z) has the two vectors ^-Bj- j . 

(-^i?) linearly independent in . {Hence, given the factorization 
m'f^^ (z) — zm\^^ [z), cyclicity of the reduced subspace 
({-2, -3,..., -2.9 + 1}), 

i.e., 

spanjz"'' ; 2 < k < 2g - 1} , 
holds for a generic subfamily {B} in 7'g(T, U2(C)).) 
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l(0) 



Proof, (a), (i) ^ (ii): If Xo{A) = 0, then Al^ = 0, and therefore 

(7 TT\ _ 4(0) 4(9) 4(0) 



i.e., A. 



ia) 
1,1 



0. This means that the coefficient matrices in the expansion 



A{z) = + ^(I'z + • • • + A'^^ha 
satisfy the conditions in Figure and, if we define matrices 



by ( |7.5| ) above, then it follows that A{z) — zB{z) where 
(7.12) = + B 



Hence (ii) holds. 

(ii) (i): This is clear from reading (7.£) in reverse. 

The equivalence (i) <^ (iii) follows from the observation that the following sum 
representation 



eo 



holds for some n and G ({—1, —2, . . . }) if and only if 

K---KK^o^{{~l,-2,...}) foralHi,...,i„. 
The conclusion can therefore be read off from the following general fact: 



It • • • Tt^t eo e <!o ■ • ■ At]o^o + ({-1, -2, . . . }) . 

(b), (i) ^ (ii): If Xo{A) = 1, then A^g^ = for fc > 0, and conversely. This 
follows from the identity 



(7.13) 



^(0)*^(0) _| ^ 



1 
1 



which is part of the defining axiom system for A. Hence, the result follows from 
|BrJoOO , Theorem 6.2], once we note that the only polynomials b{z) such that 
\b{z)\ = 1 for all z S T are the monomials; sec also [BrJoOO, Lemma 3.1]. 

(c): We already showed in Section ^that (i) implies the first of the conditions in 
(ii). The second one then follows from (7.11), i.e., the second line in (7.10) follows 
from the first one. 

(ii) (i): This follows from (b) above. For if Xq{A) — 1, then it follows from 
(b) that 

(7.14) L\R) = [02(Ceo)] ® [O2 ({-1, 

and so eo is not in the subspace 

[02({-l,-2,...})]; 



.}>], 



and the same argument, based on (7.11), shows that e_(2g+i) is not in 

[O2 {{..., ~2g + l,~2g})], 

concluding the proof of (c). 

(d): We already saw that if Xo{A) = 0, then the conditions in (7.10) hold 
if and only if 1 is not in the spectrum of the matrix from (7.6). Having now 
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(z) — zm\ {z) from (a) above, we can use the identity (7^) relating the 



Ra{ ■ , • )-numbers to the Rb{ ■ 

if and only if 
(7.15) 



numbers. But 1 is in the spectrum of the matrix 

di c 
c do 



(i-do)(i-rfi) = icr 



The matrix on the right-hand side in (7.6) is of this form, and 



(7.16) 

by Schwarz's inequality Here 
(7.17) c:=5:^<). 



< dodi 



^0 -EK!^ 



< 1, 



^1 -EK^ 



< 1. 



But using ( [7.11D and ( |7.13| ), we also get do + di < 1. Now ( |7.15D -( [7l^ ) yield 
1 < do + di, and therefore do + di = 1. Substituting this back into formula (7.15) 
then yields dodi = |c| , which amounts to "equality" in Schwarz's inequality (7.16); 



and so the corresponding vectors (d)(iii) are proportional. We already noted the 
equivalence of (i) and (ii) in (d); and we just established that the negation of (i) 
amounts to linear dependence of the vectors in (d)(iii). So (i) is equivalent to 
the linear independence, as claimed in (d)(iii). This completes the proof of the 
proposition. □ 



Remark 7.2. Let loops A and B be as in Proposition 7T(a), see also Figure y, 
and let T^"^^ and T'^^ be the corresponding wavelet representations. Then, as a 
result of the theorems in Sections ^ and |^, we conclude that T^"^^ is irreducible if 
and only if T^^^ is. Since the factorization in Proposition 7.1(a) corresponds to 
Ao (A) = 0, we conclude that the general irreducibility question has therefore been 
reduced to the case Ao (A) > 0, which is the subject of the next section. 

8. An EXPLICIT FORMULA FOR THE MINIMAL SUBSPACE 

Given N = 2, and A G Vg (T,U2 (C)), we considered the wavelet representation 
r(^) of O2 on (T) ^f. We showed that 



(8.1) 



z 



-(2S-1) 



is TI^"^^ *-invariant and cyclic (in (T)) for the representation T'^"*^. But we also 
showed that the first one of the basis vectors, z", is in the cyclic space generated 
by 



-2(9-1) 



and the representation, if and only if = for all i. 



Specifically, setting 
(8.2) 



Ao {A) 



E 



A 



(0) 



we showed in Corollary |6.4| and Proposition |7.1| (a) that 
(8.3) eo e [O2 (e_i, e_2, • • • , e_2(g-i))] 

if and only if Ao {A) = 0. Hence, it follows that /C is not minimal (in the sense 
of the following definition) if Aq {A) — 0. In other words, if Ao {A) — JC then 
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contains a strictly smaller subspace which is both T^"*' * -invariant and cyclic. In 
this section, we show the converse implication. But first a definition: 

Definition 8.1. We say that a subspace £ C /C is minimal if it is T^"^-' *-invariant, 
cyclic for the representation T^^^ and minimal with respect to the two properties, 
i.e., it does not contain a proper subspace which is also T^"^^ *-invariant and cyclic. 

We will now prove the converse to the above-mentioned result, showing, in par- 



ticular, that if Ao {A) > 0, then /C is generically minimal; see Corollary 8.7. 

Theorem 8.2. Let A £ Vg (T, U2 (C)), and let T^^^ be its wavelet representation. 
Let K, — {z^, z^^, . . . , z~^'^^~^'^) as in ( |8.l| ), and assume Ao (A) > 0. Then K, con- 
tains a unique minimal subspace C, i.e., C is T'^^^ * -invariant, cyclic, and minimal. 
Lt is spanned by the complex conjugates of the following family of Ag functions: 

(8.4) Aij (z) z''^^ , where i,j € {0,1} and k e {0,1, g - 1} . 
(Here k varies independently of both i and j .) Moreover, within the class 

(8.5) Vg{T,\J2{C)), Xo{A)>0, 

the dimension of C is 2g, for a generic subfamily, and so, for this subfamily, C = K,, 
and IC itself is minimal. 

Remark 8.3. An immediate consequence of the definition of the subspace 
£ C /C C L^(T) is the following formula for the "deficiency space": 

(8.6) /Ce/: = /\ker(PT,P) 



where P denotes the projection onto /C. We will show below, using (B.6) and 



Corollary 5.3, that £ = /C if Ao (A) > 0. This means that IC itself is then the 



unique minimal subspace when the loop A does not have m["^^ (z) e zC [z], as in 
Section ^. But before arriving at the conclusion, we must first derive several a 
priori properties of C 



Proof of Theorem 8.2. The details are somewhat technical, and it seems more prac- 
tical to first do them for the special case when g = 3, and then comment at the end 
on the (relatively minor) modifications needed in the proof for the case when g is 
arbitrary g > 2. 



Using the terminology of {5A), we then get 

(8.7) ^(z) = + .1(2)^2^ 

where A^"^^ ^ and A*^"', A^-^\ and A^^^ are 2-by-2 complex matrices satisfying 

2 

(8.8) =Vil2. 

fe=o 

When A is given, we denote that corresponding wavelet representation by T^^\ or 
just T for simplicity. Recall 

(8.9) (TJ)(z) = Y,A,,{z^)z^f{z'), 

j 

or simply 

(8.10) T,/(z) = m(^)(z)/(z2), 
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where 
(8.11) 



m 



(A) 



(^) 



7 . J' 



(z^) Z^. 



As we saw in (p.l[), the subspace 

(8.12) JC = (z", 0-3, z--*, 

is then r*-invariant, and also cyclic for the representation. But the issue is when K, 
is minimal with respect to these two properties. The minimality of some subspace 
C G K, then means that C is T*-invariant, and cyclic, and that no proper T*- 
invariant subspace of C is cyclic. 

In working out details on £, we use (S.9)-(3.11) in conjunction with (5.1C)-(5.11), 
and it is more helpful to work with the complex conjugates 



(8.13) 



M:=C = 



and so M. consists of polynomials of degree at most 5. It follows from (|8.9|)-( pjT ) 
that M is then spanned by the functions (polynomials) in the following list: 



^.14) 



A, J (z) z 



fc+j 



i = 0,l, J -0,1, fc-0,1,2. 



Hence, by ( S.13 ), L consists of the space spanned by the complex conjugates of 
these functions. By (^.12|) it is clear that £ C /C. 



The proof of Theorem B.2 will now be split up into several lemmas: 
Lemma 8.4. The space C is T* -invariant. 



Proof. Now, the functions AiAz) in ( S.14 ) are the matrix elements of the loop 
A(z), and so it follows from (8/7) that each of them is a polynomial of degree at 
most 2, say 

(8.15) a (z) = Co + ciz + C2Z^ 
(since A (z) itself has degree 2 when g ~ i). Hence, 

(8.16) T* (a) = coA,,o (z) + ciA,,i {z)z-^ + C2A,,o (z)z-\ 
or equivalently, 
(8.17) 

Now set h (z) 



), we then get 



T* (a) = coA,,o [z] + ciA,,i (z) z^ + caA^^o {z) z^ 
z'^Pa (z) where a is as in ( ^.15| ). From ( 
T* (b) = z-'^T* (a) , 

or equivalently, 

(8.18) T* (b) = zP (coA,,o (z) + ciA,,i (z) z + c2A,,o (z) z) , 

using (B.17). In view of (8.14), we need then only to compute the following: 

T* [za (z)^ = coAi,i {z)z~^ + ciAj,o {z)z~^ + C2Ai,i (z)z-^, 
or equivalently. 



[za (z)^ = coAi^ (z) z + ciAi^o {z) z + C2Ai4 (z) z^ 



Now putting the formulas together, we get the value of T* on each of the functions 
( ^.14 ) which go into the definition of £, and the conclusion of the lemma follows. □ 
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Lemma 8.5. The space C is cyclic. 



Proof. Since both £ and /C are r*-invariant, the conclusion wiU foUow if we check 
the inclusion 



(8.19) 



/C C span (TiC) . 



For the space on the right-hand side in ( |8.19| ), we shall use the terminology [Oj/^j , 
and similarly, the space spanned by all the spaces 



.20) 



Til 



•T £ 



will be denoted [O2C] . In ( B.20 ), we vary the multi-index (ii, 12, . . . , in) over a ll th e 
2" possibihties. It follows from the T*-invariance of JC (in (^)) and £ (in (^.4[)) 
that we get different families of nested finite-dimensional subspaces: 

(8.21) K. c [OllC] c [OlK] c [OllC] c • • • c [O^K.] c [O^+^IC] c • • • , 
and a similar sequence for £. Since /C is cyclic, we have 

(8.22) y [O^IC] = (T) (= f). 

n 

But K. C [O2C], for some n, so £ is also cyclic. The conclusion of the lemma follows 
from 



.23) 



.21) and (8.22), once we check that 



and so n = 1 works, and[C'2/C] C 



for all p. 



Turning now to the details: Since /C is spanned by z fc = 0, 1, . 
check that each of these basis functions has the representation 

(8.24) z-'^^ym 



. , 5, we must 



for loji G £, where we refer to ( B.14| ) (see also ( ^.4[ )) for the characterization of 
the spa ce £, or rather A4 := £. 
But 



(8.25) 



.24 ) is equivalent to the assertion that 



for all i — 0,1 and fc = 0, 1, . . . , 5; and ( B.25 ) can be checked by a direct calculation 



which is very similar to the one going into the proof of Lemma 3.4. Specifically, 
using ( ^.1C| )-( 5.11 ) we get the following: 

T* (z") = e £, 

= A, I (z)z-^ e £, 



(8.26) 
(8.27) 
(8.28) 
(8.29) 
(8.30) 
and finally 
(8.31) 



t: 

t: (z-) 

T* iz-^) 



= A 



A^^n{z)z 1 e £, 
A~^z-^ e £, 



A^,l iz)z-' e £. 



Recall that the complex conjugates of the functions on the right-hand side in this 
list are precisely the ones from ( 8.14 ), or equivalently, (^.4[). This proves ( |8.23 ), 
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and therefore the cychcity of £, which was claimed in the lemma. As a bonus, we 
get from ( ^.26D -( |8!3ll ) that the inclusion £ C ({-1, -2, -3, -4}) holds if and only 
if \q{A) — 0. To see this, use the fact (for 5 = 3) that 



A 



(2) 



□ 



Lemma 8.6. The space C is minimal in the sense of Definition 



.1 



Proof. We will establish the conclusion by proving that if £i is any subspace of K, 
which is both r*-invariant and cyclic, then £ C £i. So in particular, C does not 
contain a proper subspace which is both r*-invariant and cyclic. 

Now suppose that some space Ci has the stated properties. Since it is cyclic, we 
must have 



.32) 



satisfied for some n £ N. As noted in the proof of Lemma 8.5, this is equivalent to 
(8.33) T*^---T*T*^{z-'')e£i 

for all ii, . . . ,in G {0, 1}, and all k g {0, 1, . . . , 5}. But we also saw in the proof 



of Lemma B.4 that the functions on the left-hand side in^ 
which are listed in ( ^.14 ). Note that the functions in (8.1 
given by 



33| ) are precisely those 
or (8.4), are those 



(8.34) 



1 = 0,1, 



A: = 0,1,. 



But Ao {A) > by assumption, so for some i, we have Af^ ^ 0, and the calculation 



in the proof of Lemma 8.4, and in the previous two sections, then shows that 
the families of functions in (B.34) and ( 8.33| ) are the same, i.e., we get the same 
functions in (8.33) for n > 1 as the ones which are already obtained for n = 1 in 
( ^.34 ). This is the step which uses the assumption Ao {A) > 0. Since C is spanned 



by the vectors in (B.34 



worked out in Remark B.8 below 



), the desired inclusion £ C £i follows. 



More details are 
□ 



Proof of Theorem 8.2 concluded. The result in the theorem is now immediate from 



the three lemmas, and we need only comment on the size of the genus g. We argued 
the case g = 3; but, for the general case, IC is spanned by z~'^, k — 0,1, . . . ,2g — 1, 
and the functions from the list 



.35) 



.14), or equivalently 
A, (z)z'=+^ i = 0,l, j = 0,l, fc = 



0,1 



4), will then be 
1. 



,g- 



Oth erwise, all the arguments from the proofs of the lemmas carry over. See Remark 
B.8 for more details. □ 



(0) 



i?(0,0)oo > 0, then 



Corollary 8.7. When g is given, and Ao := J2i ^ 

C = K for a generic set of loops A in Vg (T, U2 (C)). 

Proof. The proof comes down to a dimension count. Since K, — (z", . . . , 2^(29-1)^ 
is of dimensi on 2g , we just need to check that the space C (c IC), spanned by the Ag 
functions in (8.35), is of dimension 2g for a generic set of loops A in Vg (T, U2 (C)), 
and that can be checked by a determinant argument based on the conditions for 
the matrices A^^) 



,A^s-i) defining A (z); see (0)-(ra above 
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The above-mentioned dimension count is based on the following consideration 
(which we only sketch in rough outline). A possible linear relation among the 
functions from (S.4) takes the form 

(8.36) 5]^5]a,,,fcA,,,(z)z^+'=^0, 

i j k 

where the i, j summation indices are 0, 1, and the k summation is over 0, 1, ... , <?— 1. 
As a result, we get the following system of relations: 

1 



.37) 



Ci^j^kAij (z) = (mod; 



i=0 



for all j = 0, 1, and all fc = 0, 1, . . . , 5— 1. Note that (8.37) is a matrix multiplication. 
Using finally 

■ i(0) ^ 



(8.38) 



>0, 



we see that the dimension of the space spanned by {Aij {z) 2:^+'''} is 2g, as claimed. 
See Remark B.8 and Observation fe.9| for details. □ 



Remark 8.8. A more detailed study of the space C will be postponed to a later 
paper, but one point is included here: The function z'^ (= cq = 1) is in C if and 
only if the polynomials 



(8.39) 



{^,0 {z),Aj,j{z)z}^ 



do not have a common divisor. This follows from ( 8.26| )-( ^.31 ). Indeed, for Cq to 
be in £, we must have the existence of hij (z) € C [z] such that 



(8.40) 



But by algebra, this amounts to the assertion that the family of polynomials listed 
in (B.39) is mutually prime within the ring C [z]. Also note that, by the result in 
Section 7, monomials such as d (z) = z are not common divisors in the polynomials 
of (B.39) if do (A) > 0. In fact, a possible common divisor d{z) ^ C [z] for (B.39) 



yields the following factorization: 
(8.41) A,^o{z) = d{z)h^o{z), 

{ki^o (z) , kj,i (z) £ C [z]). Hence: 



Aj.i (z) z = d (z) fcj.i (z) 



Observation 8.9. If do (A) > 0, then eo € C. 



Proof. For if not, the greatest common divisor d (z) of the family (B.39) would have 
a root 7 e C \ {0}, i.e., d (7) 



0. By 



..42) 



^(7) 



41), we would then have 
0, 



where A eV (T,U2 (C)) is the originally given loop. Recall that, by (p.lD, we may 
view A{z) as an entire analytic matrix function, i.e., an entire analytic function, 
C M2 (C), whose restriction to T takes values in U2 (C). (These are also called 
inner matrix functions [ PrSe86| .) But (B.42) is impossible (for 7 7^ 0) in view of 
Lemma 5.1 and its corollary. We will give the details ioi g — 3, but they apply 
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with the obvious modifications to the general case of 5 > 2. If dsll ) holds, then 
by Lemma 5.1, 

(8.43) V-^A (7) = {Q^ + 7Q0) [Qi + iQi) = 0. 

where we use the projections Qq, Qi in from (^). Setting (z) = Qj~ + zQj, 
j = 0, 1, ( 8.43 ) then yields the following estimate: 

(8.44) = A, (7)* Ao (7)* Ao (7) A, (7) - (7)* (q^ + Qo) A, (7) 

> min (1, |7|') • {A, (7)* A, (7)) = min (l, I7I') • (q^ + I7I' Qi] 

> ( min ( 1, I7 



•I2, 

where the order > is that of positive operators on C'^. But this is impossible, since 
7 7^ 0. The latter is from the assumption do (A) > 0. □ 



This last argument in this proof also serves as a proof of Corollary 5.3, and this 
corollary is again the basis for the following stronger result, which we now sketch: 

Theorem 8.10. If Xq (A) > 0, it follows that C = IC. 

Proof. By Corollary and ( |3.12 ), we have the estimate 

(8.45) ^|mi(z)|^ > (min (1,121^))" (l + for aU z G C. 

i=0 

If fc e /C e £, then by (|!|), we get 

(8.46) Pm, (z) A: (z^) = 0, i = 0, 1, z e T. 

If Ao {A) > 0, then the value z = is not a common root of the two complex 



polynomials mo (z), mi (z), and so by (8.45) the two polynomials mo and mi have 
no common roots at all, by Proposition 7.1(a). Therefore, when (8.45) and (3.46) 



are combined, we get fc = 0. Hence K. Q C — 0, and the proof is completed. 



□ 



Remark 8.11. Not everything that works out easily in the case g = 3 generalizes 
immediately to g > 3: The case g = 3, and arbitrary N , amounts to a choice of 
two projections, say P and Q, in C^, and the finite-dimensional representations 
of the algebra generated by two projections are completely known by folklore (see. 



e.g., |JSW95]). In fact, it can be easily checked that this is the same as displaying 



the finite-dimensional representations of the Clifford algebra with two generators, 
A, B, say. The relations between A, B are: 



(8.47) A* = A, 
For any such pair, set 

(8.48) P = 



B* = P, 



A' 



B^ = 1 



N, 



and AB + BA = 0. 



1 



1 + A~B), 



Q 



1 



(t~A 



Then it is immediate that P and Q are projections, i.e 
Conversely, if P, Q are any projections, set 

(8.49) A^P-Q, B - P -Q, 

and an easy calculation shows that 



P), 
P 



etc. 



47) is then satisfied. Since the finite-di- 



mensional representations of (B.47), the Chfford algebra C2, are known |JSW95|, 
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we then get a useful classification of 1^3 (T, Ujv (C)). But these comments do not 
carry over to the case g > 3. (A good reference on the Clifford algebra and its 
representations is [LaMi89, Ch. 1, §5].) 

It is perhaps a little early to identify regions in the parameters 0, p where the 



scaling function x 



(x) is regular, and where it is not, but a primitive test 



would be the vanishing-moment condition of Daubcchics | Dau92 , Ch. 6] . We would 
look for values 9, p such that m'^'''^ (z) is divisible by (1 + z)^ and by (1 + . 
These are the conditions which ensure that 

fe 



(8.50) 

starting with 

(8.51) 

(8.52) 



^(Ok=o = 0, fc = 0,l. 



= (0) 



^ {x) dx, 



= / x'^i) {x) dx. 



Here ( p.51| ) is automatic since mg (—1) = 0. Recall the coordinates z 

k 



^ G K. The second one (B.52) corresponds to ( ^ j ip {^) \^=o ~ 0, or alternatively, 



klLzl (■^) U=-i — 0' ™ another form, the condition that (1 + z) 
factor of mo (z), etc. 



fe+i 



is a 



Proposition 8.12. (a) The polynomial 1711^ (z) is divisible by (1 + z)^ (see ( 8.52 )) 
if and only if 



.53) 



cos (29) + cos {2p) = 



(shown as curves in the four corners of Figure ^ in the Appendix). 
(b) The polynomial ni"^'''^ (z) is divisible by {\ + z)^ if and only if 



(8.54) 
i.e., when 

9 = cos' 



cos 29 + cos 2p — 



1 



%d sin 29 + sin 2p = 2 sin [29 - 2p) , 



-1 4 



0.89 w 0.28 TT, 



p = cos 




0.39 « 0.12 TT 



or when {9,p) is related to this pair by {9,p) —> {9 + rmr, p + nn) , m,n £ or 
{9,p) {—9,—p), or both. {The points in [0, tt] x [0,7r] are shown in Figure ^) 



_ / ma(z) 



Proof, (a) In this example, N ~ 2 and g = i . Let m (z) 
QMF-polynomials in z, viewed as a column vector. By Lemma I 

(8.55) m{z)^V [P^ + z^P) {Q^ + z^Q) a (z) , 



be the usual 
we have 



where V" £ U2 (C), P and Q are projections, and as before , a (z ) :— (l). Here 



and P, Q are the 0, p projections specified in ( 6.50 ) (see also ( A. 11 ) 
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below). We need only check that j^m (z) \z=~i = ( ) for some number ^q. But 
(8.56) (z) = V2zP {Q^ + z^Q) a{z) + V {P^ + z^P) 2zQa (z) 



V{P^ + z'P){Q^+z'Q) [ 5 



Substitution of z = — 1 yields 

(8.57) m'{-l) = -2V{P + Q)(^ \^+v(^l 
Hence 

(8.58) ^(^ + ^)(-l)-(?)^(4 
for some number ^j^, and therefore, with 

(f>.^Q\ P n-f ^ -f cos 26* + cos 2p sin 26* + sin 2p 

(8.59) P + Q-]^^ ^ j + 2 sin26' + sin2p - (cos 26* + cos 2p) 

we have 

2 \ / cos 26* + cos 2p sin 26* + sin 2p \ / 1 \ f \ _ f 



2 y \^ sin 26* + sin 2p - (cos 26* + cos 2/?) J \ -'i^ J \ J \ ^Ci 

The result (a) follows. 

Part (b) follows upon solving 

(8.60) ^miz)U^^^=(^ 

for some number ^2 1 in addition to the conditions in (a) . But we only need to work 
out the next derivative, using ( 8.53| ) to eliminate the cosine terms where possible: 

(8.61) y-i^m(z)|,_i 

= (2 (P + Q) + 8PQ) \ ^-4(P + Q)(^ 5 

|-3(sin20 + sin2p) \ / 3 - 2 (sin 26* + sin 2p) 
-| + (sin26l + sin2p) )^\ -1 + 2 (sin 26* + sin 2/9) 



+ 2 



cos2 {9 - p) + sin2 {e - p) 
sin 2 (6* - p) - cos 2 (6* - p) 

0]^( ^3 
J \ -^3 

where we used 

^ff ^ ^\ / cos26' + cos2p sin 26* + sin 2/9 



^'^^ 4 VV 1 y ' V sin 261 + sin 2p - (cos 261 + cos 2/9) 

cos 2(9- p) - sin 2{9 - p) 
sin 2 (6* - p) cos 2 (6* - p) 

Noting again that the right-hand side of ( |8.61 ) takes the form (^-5)1 arrive at 

(8.63) -2 (sin 261 + sin 2p) + 4 sin 2 (g - p) = 0, 
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which directly gives the second part of (S.54) 



The specific sohition is found by transforming (8.63) into 

(8.64) (1-2 cos 2p) sin 26* = - (1 + 2 cos 29) sin 2p. 
Squaring and using ( ^.53 ) to ehminate p yields 

(8.65) 2 cos^ 20 + 4 cos 26* 
which gives 



3 
4 



0, 



(8.66) 



cos 29 



by the quadratic formula, using |cos2(?| < 1 to choose the positive radical. Working 
this back through (B.53) yields 



(8.67) 



and substitution of 



cos 2p = - — 



6(;) and (8.67) into (8.64) then shows that sin 26* and sin2p 



must have the same sign, so that the solutions stated in the proposition are the 
only ones possible, the numerically exhibited pair being those for which sin 20 and 
sin2p are both positive. □ 



Remark 8.13. The Holder-Sobolev exponent is at least as good as 0.84 when we 
have a (1 + z)"^ factor of ■mQ{z), and at least as good as 1.136 if (1 + z)^ is a factor 
|LaSuOO| [|Vol95| ; see Figure H. 



Appendix (by Brian Treadway) 

In this appendix we show some cascade approximations of wavelet scaling func- 
tions for the example with g = 3 discussed in Sections || a nd |^ abo ve. The local 
method of cascade iteration works here just as it did in prJo99bt . It is just a 
matter of enumerating terms. 



The direct method of iterating the relation (1.7), i.e., 



2s- 1 



(A.l) 



(p{x) = ^ akf (2x - k) 



k=0 



proceeds by translating a distance k to the right, multiplying by a^, summing over 
fc, and scaling down by 2. This takes an expression in which every term has n 
factors, each of which is an ak for some k, and turns it into an expression in which 
every term has n + 1 such factors. In fact, every ordered product of n a^'s occurs 
exactly once, at some dyadic point in the n'th stage. Multiplication of these factors 
is commutative, but we forget that for the moment and take care to add the new 
factors at the left. 

For example: With 5 = 3 there are 6 coefficients, oq, ai, . . . , as. 

n = 0: 1 point, 1 = 6" term, no factors of a^, so the vacuous term is just 1, and 
we have the Haar function: 



a; = Lp = 1 
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Figure 2. The ultra-smooth wavelet scaling function: Lpg ^ [x) for 
9 « 0.284 TT , p K , 0.124 TT, with two vanishing moments of "0 (see 
Proposition 8.12| (b) and the discussion preceding it). The Holder- 
Sobolev exponent of this <p is > 1.136; see Remark |8.13 . 



1: 6 points, 6 = 6^ terms, as follows: 







X = 


(fi = ao 


X - 


= 3/2 




= 03 










X = 


1/2 if = ai 


X - 


-- 2 




— 04 










X = 


1 (fi = 0.2 


X - 


= 5/2 




= 05 






n — 


2: 16 points, 36 = 


& terms, as follows: 












X 


= 


if = aoao 




X — 


2 




— 0400 H 


^ 03O2 - 


V 02O4 


X 


= 1/4 


if = aofli 




X = 


9/4 




= 04O1 - 


^ 03 03 - 


1- 02O5 


X 


= 1/2 


ip = amo 4 


0002 


X = 


5/2 


f 


= 0500 H 


^ 0402 - 


f- 0304 


X 


= 3/4 


(p — aifli 4 


oooa 


X = 


11/4 




= 0501 H 


- 04O3 - 




X 


= 1 


f = 0200 + 


aia2 + 0004 


X = 


3 




= 0502 H 


- 04 04 




X 


= 5/4 


if = 0201 4 


oioa + 0005 


X ~ 


13/4 




= 0503 H 


- 0405 




X 


= 3/2 


if = aaflo 4 


0202 + O1O4 


X — 


7/2 




= 05O4 






X 


= 7/4 




02O3 + 0105 


X = 


15/4 




= 0505 






n = 


3: 36 points, 216 = 


: (r" terms, grouping 


left as 


an 


exercise 


for the reader 



One could write down all 6" terms at the outset of each stage, and then just ask 
which ones go with which of the values of x. The answer is: 

(A.2) 

The indices of the o's, interpreted as digits at the right of the fraction point 
in a non-unique positional binary number system with six digits instead of 
two, give the x to which a given term should be assigned. 
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As an example from the list above, the term 0401 goes with x — 4 (2~^) + 1 (2^^) — 
9/4. All the others (at any cascade stage) can be assigned in the same way. 

In the ordering used above, if we go from one stage to the next by the "translate 
by k, multiply by a^, sum over fc, and scale down by 2" method, terms will be built 
up by adding factors from the left, while if we use the "local linear combination" 
method from | BrJo99b| Appendix] and |Dau92, Section 6.5, pp. 204-206], terms 
will be built up by adding factors from the right. In either method, the full set of 
6" terms (or 4" in |Br Jo99l| , where g — 2) will be obtained, without duplication, 
and each term will be assigned to the same x regardless of which method is used. 

At cascade stage n there are values of (p assigned to values Xi of x that are 
consecutive integer multiples of 2~" ranging from to (g'"' — l) -2"". The iteration 
begins with the ordered list of values of (p at those points on the x-axis: 



(A.3) (^(")(0),/")( 



1 • 2-") , vjf") 



2 • 2- 



((,<-. -1) 



• 2- 



Sets of g adjacent values in this list go into the computation of each point at the 
next cascade stage, which is on a finer grid of g("+^) cons ecut ive integer multiples of 
2-("+i). Each set oi g adjacent values of ip from the list ( |A.3| ) yields 2 values on the 
2-("+i) grid by a linear combination with two sets of alternate a^'s as coefficients. 
This can be expressed as a matrix product (shown here in the case g = 3): 



( V 



(n+l) 



(A.4) 



(0) 



^(„+i) ^2 . 2-("+i)) 
(^("+1) (4 • 2-("+i)) 
<p("+i) (6 . 2-("+i)) 



^(n+i) (3.2-(«+i)) 
(5.2-("+i)) 
(7.2-(«+i)) 



\ 



^(„+i) ((g(n+i) „ 6) • 2-("+i)) ((g("+i) - 5) • 2-("+i)) 

^(„+i) ((g(n+i) _ 4) . 2-(«+i)) ((g(«+i) - 3) • 2-("+i)) 

V ((g("+i' - 2) • 2-("+i)) ((g("+i' - 1) • 2-("+i)) J 















(^(") 


(0) 


\ 









(^(«) 


(0) 




(1 • 2-") 




(0) 






(^(«) 


(1 • 2-") 




(2 • 2-") 




^(") (1 • 2- 








(2-2-") 




(3 • 2-") 




^(") ((g(«) 


-3) 


• 2-") 




((qt") - 2) • 2-") 




((g(«) - 1) 




^(") ((g(«) 


-2) 






((qf") - 1) • 2-") 









^(") ((g(») 


-1) 


• 2-") 












/ 



Thus a (g^"^ + (5 — 1)) x g matrix (the first factor on the right in (A^) above) 




partitioned from the list (A.3) of g^^"^ points at stage n, yields a (g*^"^ + (g — 1)) x 2 
matrix (the left-hand side of (A.4)), which is then flattened out to give the list of 
values for stage n+l. The number of points in the new list is the total 



number of entries in that (g^"^ + (g — 1)) x 2 matrix. 



(A.5) 



g("+i) - 2 (g(") + (.g - 1)) 
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and with g'"-* — 1, this recursively yields the number g^") of points at each stage as 



(A.6) 

In the case 5 = 3, this is 
(A.7) 



(2ff-l).2"-2(5-l). 
= 5 • 2" - 4. 



The three steps in the local method, (1) partitioning a list into rows of g ~ 
3 points (first adding 5 — I zeroes at each end), (2) matrix multiplication, and 
(3) flattening the resu lting matrix back into a single list, are easily implemented in 
Mathematica |Wol96| . AU that then remains to compute a cascade approximation 
of a wavelet scaling function ip is to specify the numerical values of the coefficients 
Uk and repeat the procedure n times, starting with the one-element list 



(A. 



(!)• 



The same local cascade relation expressed in (A. 4) as giving two values of (^1^"+^^ 
from three values of (^'^"^ can also be set up to give four values of (p^"^^'> from 
four values of (p'"-*, or five values of (^5^"+-'^) from five values of ip^"\ by combining 
overlapping ranges of the initial and final lists. The terms involved on the successive 
cc-grids are shown in the diagrams below, and the corresponding matrices are given. 









04 


as 










( 


a2 


03 


) 


(as 


in (/ 




ao 


ai 












04 


a, 5 














02 


03 


04 


as 








ao 




02 


03 






V 








ao 


Ol 


) 






05 


















0.3 


04 


05 












ai 


0-2 


03 


04 




05 







ao 


Ol 


02 




03 


V 











ao 




Ol 


/ 


04 


05 















a2 


03 


ai 


05 









ao 


Ol 


02 


03 




04 










ao 


Ol 




02 


V 
















oo 



All of these represent the same calculation of the {n + l)'st cascade stage fron 
the n'th stage: they merely collect different locally related sets of points in the 
successive stages. The advantage of the matrices that are square is that they allow 
successive stages to be expressed as powers of the matrix. 

An eigenvector decomposition yielding an explicit expression for the n 00 limit 
of the cascade stages, like that done in | BrJo99b| , Appendix], could be done here 
using the 4x4 matrix above (here the matrices are written to act on the left). The 
detailed calculation is much more extensive than in the 2x2 case of iBrJo99b| , and 
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we will not present it here, but note only that the two 5x5 matrices above have sim- 
ple eigenvalues and left eigenvectors (in addition to the row ( 1 1 ... 1 ) that 
all the matrices have): for the first, eigenvalue 05, eigenvector ( 1 ), 
and for the second, eigenvalue oq, eigenvector ( 1 ). Since the start- 
ing list for the cascade computation is ( 1 ), to be "padded" with zeroes on the 
left and right, these two eigenvectors occur explicitly at the first cascade stage and 
continue thereafter. This shows that 1^9^"^ diverges at one end of its support interval 
or the other when one of these eigenvalues is greater than 1, growing like when 
as > 1 or like aJJ when oq > 1. The full eigenvector decomposition then has a 
term, generically nonzero, that grows in the same way for other points Xi, and so 
the same divergence occurs at points throughout the support interval. The regions 
where oq or 05 is greater than 1, leading to this divergence in the cascade iteration, 
are shown by shading in Figure ^. 



The coefficients here are not restricted, as they are in | WanOO | , to < 1, so we 
do not have cycles in the cascade iteration: the terms grow indefinitely. As a result, 
the condition as > 1 or oq > 1 is sufficient, but not necessary, for divergence. 

For the example with g = 3 specified in (|6.50| ) above, the formulas for the 
coefficients ao,ai, . . . , 05 of the polynomial m^'^ (z) in (B.ll) may be derived as 
follows. From ( 6.43| ) we have 



(A.9) 
with 

(A.IO) 



A{z)^V{Qi+zQo){Qj + zQ) 



V = 



1 1 

1 -1 



(A.ll) 



cos U 
COS 9 sin ( 



cos f sm ( 



and 
(A.12) 

Then the coefficients aq, ai. 



1 
1 



Qe' — Qe+{TT/2) 
, 05 are: 



cos 20 sin 29 
sin 20 - cos 20 



(A.13) 

ao = 




cos 20 - 


- sin 20 


— cos 2p 


— sin 2p - 


- cos{20 - 


2p) + 


sin(26l - 


2p)), 


ai = 




cos 20 - 


- sin 20 


+ cos 2p 


— sin 2p - 


- cos(20 - 


2p)- 


sin(26l - 


2p)), 


a2 = 




cos(26' 


-2p)- 


sin(26l - 


2p)), 










0.3 = 




cos(26' 


-2p) + 


sin(26l - 


2p)), 










04 = 




cos 20 - 


f sin 20 


+ cos 2p 


+ sin 2p H 


- cos(20 - 


2p) + 


sin(26l - 


2p)), 


05 = 




cos 20- 


f sin 20 


— cos 2p 


+ sin 2p - 


- cos(26' - 


2p)- 


sin(26l - 


2p)). 
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These can be seen to meet the conditions (1.5)-(1.6) for the coefBcients of a scahng 
function, The even- and odd-indexed coefficients also sum to a constant separately: 

2 2 

(A.14) ^a2i = ^a2»+i = 1 

i=0 i=0 



(see I ReWe98, eq. (9.12)]), which is what makes the constant vector ( 1 1 ... 1 ) 
an eigenvector or the a^-matrices above. Of course, the coefficients a^, as functions 
of 9 and p, have the periodicity (with period tt in both angles 9 and p) of the 
projections Qe and Qp they were derived from: 

(A. 15) ai {9, p) = ai {9 + rmr, p + mr) , m,n€'Z. 

In addition, they are related in pairs by the reflection relation 

(A.16) a,{9,p) = a5^,{-9,-p), i^0,...,5. 

These relations carry through the successive stages of the local cascade computation 
(A.4) as the periodicity 

(A.17) ^S(^)=^1"+Up+«.(^) 
and the reflection symmetry 

(A.18) (x) = fi%^^_^ (.) = (4") - .) , 

where Xf is the last point to which a value is assigned by the ?ith stage. By (A. 7) 
we have 

(A.19) x^"^ = fq(") - l) • 2-" = 5 - 5 • 2"". 



f 

There are pairwise relations involving a translation in the (0, p) plane by half the 
period: 

(A.20) a2.+,(0,p)-a2(2-»)+j(f?-|,p-|), « = 0,1,2, j=0,l. 

There are also twofold and threefold affine symmetries, such as the invariance of 02 
under the twofold transformation 

(A.21) 9^-9, p^-p+1^ 

or the invariance of oq under the threefold transformation 
(A.22) e^-p+^, p^9-p+^, 

which has ap's three local extrema as its fixed points in the 7r-periodic context. 

Contour plots of oq and a2 in the 9, p-plane are shown i n Fig ure ^ below; the 
other fli's can be de rived from these by the translation in ( |A.20[) or the rotation 
around the origin in ( A.16 ), or around the point (§, f-) when this is combined with 
translation by tt in both angles. 

Discussion of what variation has between Xj and x^+i is somewhat meta- 
physical: it only matters that <y9*^"' (xi) is associated with the point Xi = i ■ 2^". 
For example, if we said (^(") (x) = v'^") {i ■ 2"") for i ■ 2-" < x < {i + I) ■ 2^" and 
{x) = for (« + 5) ■ 2-" < X < + ■ 2-", then the same "shape" would 
occur in every bottom-level interval of every cascade stage, preventing continuity 
from appearing in the n ^ 00 limit even if it would otherwise have appeared. Other 
variations within the shortest dyadic intervals at a given stage have a similar arbi- 
trariness; iy9^"^ is really only defined pointwise for finite n. When n goes to infinity. 





Figure 3. Contour plots of ao and 02- For uq, the three hnes 
= 7r/4, p — (or tt), and 9 — p — 7r/2 are contours of oq = 0, 
whose intersections are saddle points; all the maxima and minima 
of ao are located on the dashed line = 7r/4 (or 57r/4), displayed 
in the lower plot. For 02, contours of constant 02 are diagonal lines 
of the form 9 — p = constant. 



continuity arises in some cases. On the other hand, for a given x and finite n, 
(x) as a function of 9 and p is continuous, since it is a polynomial (of order n) 
in the a^'s; when the order n goes to infinity, singularities arise. 

To decide when the corresponding scaling function ^pe,p (x) generates a wavelet 
in the strict sense or merely a tight fra me, as dis cussed i n Section |l] above, we use 
the method of Cohen [|Coh92b| , |CoRy95|| (see also jBEJOOl ): We identify cycles on T 



for the doubling map z 1— > z^, i.e., a finite cyclic subset unequal to {1} and invariant 
under z 1-^ z'^. The result is that (p (x) generates a "strict" wavelet if and only if 



(A.23) 



zeT 



{-z) = 0} 



does not contain a nontrivial cycle. 
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The cycles on T are not subgroups of T but rather cycUc orbits on T under the 
z 1-^ z"^ action of one of the cychc groups Z^, fc = 1, 2, . . . . Such a cyclic orbit Ck 
with k distinct points zi , . . . , Zfe must be of the form zi — > Z2 ■ • • — > 2^ — > zi , 
where z^_^_i — zf if i = 1, 2, . . . , fc — 1, and z^ = z-^. Hence points c in an orbit 
Ck must satisfy c^*" — c, and each c must be a (2*^ — l)'th root of 1. Different 
orbits must be disjoint, and their union will be invariant under z i-^- z^ acting on 
T. The converse is not true. Note also that we can have different (2*^ — l)'th roots 
c of 1 defining different cyclic orbits for the same k. Iffc=lor/c = 2, then 
in each case there is only one orbit, but if /c = 3, there are two choices. Since 
TOq^''''* (z) is fo r each 0, p a, polynomial of degree at most 5, the cardinality of a cycle 
contained in ( A.23| ) is at most 4. Thus, if z is contained in such a cycle, we must 



have one of the possibilities z'^ — z, z"^ — z, z^ — z. Hence the cycles of length 
at most 3 are {1}, {uj,uj'^} where u = e*^'^/^, {C,C^,C''} where C := e'^'^/'^, and 
{C, C^,C'*} — {C^iC^iC^}- But as r7io(— 1) — always, (z + 1) is always a factor 
of Too (z), and since the cycle should be different from the trivial cycle {1}, we are 
reduced to the case {wjW^}. The other cycles would make m^'''^ divisible by a 
polynomial of degree at least 4. 

Thus we have the following cases: tti^q'^^ (z) may be divisible by 



(A.24) p3 (z) = Y[{lu'' + z) = 1 + z^ 



A;=0 

by 

(A.25) p4 (z) = (1 + z) (C + z) (C' + z) (C* + z) 

= l + (3z- z'^ + f3z^ + z^, 

or by 

(A.26) pf'> (z) = ={l + z) {(' + z) (C^ + z) {e + z) 

= 1 + /3z - z^ + /3z3 + z^ 

where C, (as above) and (3 are defined as 

(A.27) C:=e^^, ;3 := 1 + C + + = ^ + 

Proposition A.l. There are only Jour cases of QMF-functions 

(A. 28) flQ + fliz + a2Z^ + a^z^ + a^z^^ + asz^, a^ G M, 

which give tight frames that are not strict wavelets. In addition to the V3 (T, U2 (C))- 
conditions, they satisfy 

5 

(A.29) = 2. 

k=0 

The four correspond to the three loops 

l/zzM 1/1 z\ J_/^^^l 



71 V 1 ; ' v^yz ^z^ 71 V -1 
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and the loop 

The wavelet representation T*^"*^ is irreducible for the first two of the four, and 
reducible for the last two. The values of Xo {A) are as follows: 1/2, 1/2, 0, and 1, 
respectively. The first three have cycles of order 2 and the last one a cycle of order 
4. The corresponding system of coefficients is as follows: 

ao ai a2 as 04 05 

~o i i~ ] 

(A.32) 1 1 \ two-cycle 

10 10 J 

1 1 four-cycle 

and so all four cases are Haar wavelets. The scaling functions Lp [x) may be taken 
as in Table |^. 

Proof. If mo (z) is divisible by 1 + z^, then its six coefRcients oq, ai, . . . , 05 must be 
of the form cq, ci, C2, cq, ci, C2, and the associated loop T — s- U2 (C), 



(A.33) 



Cq -\- C2Z + CiZ^ Ci+CoZ + C2Z^ 

C2 + CqZ + CiZ^ ~ (Cl + C2Z + CqZ^) 



The corresponding U2 (C)-conditions then yield: 

(A. 34) 2coC2 + cico + C2C1 = 0, cqCi + C1C2 = 0. 

Substitution of the second into the first yields C0C2 — 0. Hence, of the three numbers 
Cqj Cl, C2, at most one, and therefore precisely one, can be nonzero. But each of the 
three cases is determined up to scale, and condition ( |A.29 ) decides the scale. We 



are therefore led to the three loops in (A.30), and the rule ( 8.1l|) then gives the 
three scaling functions if (x) which are listed in ( |A.32| ) and Table g. 

The cycle of the last line in Table | is of order 4. Let ^ e^^'^/s = (A 
^i2TT/i5y Then the cycle is {fi, fM^ , fi"^ , fi^} , and 0^=0 (m'' + — z^ -h 1, which is 
the mo (z) for the last line of Table ||. (It is from a root of 1 of order 2' — 1 (= 15) 
for I = 4.) 

The other length-3 loops which would be possible are, as noted, {CjC^jC''} ^i^id 
{C,C^,C'^}, with = e*^'^/''. We claim that they do not in fact occur. 

If one of them did occur, then the corresponding mo (z) would be divisible by 
either p^ (z), or by p)^' (z). But pi (1) = p^ (1) = 2, so the factorization would 
be mo (z) = Pi (z) I (z) where I (z) ~ a + (1 — a) z. (We have picked the nor- 
malization of mo (z) given by mo (1) = 2 for convenience.) From the formulas 



(A. 10) and (A. 13) we note that the coefficients ao,ai, . . . , as are real. Divisibility 
by P4 (z) means that —1, — C^, and are roots of mo (z). So the complex 
conjugates — C^, — are also roots. But that would give us all seven points, 
— C^,—C^,~C^,—C^,—C^, as distinct roots of mo (z), which is impossible 
since mo is of degree at most 5. □ 

In conclusion, when 5 = 3, the variety of the wavelets which are only tight frames 
sits on a finite subset of the full variety of all (T, U2 (C)) examples. 

On the following pages are plots, for various values of the angles 6 and p, of the 
wavelet scaling functions ipe,p {x) at the 8th cascade level, computed by the local 
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Table 2. (f{x) in the four cases: Tight frames corresponding to 
cycles of length two, and a four-cycle. 



e = 7r/4, 
P = 7r/2; 



e^O {or tt), 
p = (or tt); 



1 



1 



1 



= 37r/4, 
P = 7r/2; 



= p = w/2. 
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Table 3. Layout of scaling function plots. Round solid points (•): 
approximate locations of "ultra-smooth" wavelet scaling functions 
in relation to scaling functions plotted. Shading: divergence due 
to Co > 1 or 05 > 1. Boxes: marginal divergence due to ag = 1 or 
05 = 1. 



Il7r/12 


al 


bl 


cl 


57r/6 


ak 


bk 


ck 


37r/4 


aj 


bj 


cj 


27r/3 


ai 


bi 


ci 






p. 34 




77r/12 


ah 


bh 


ch 


tt/2 


ag 


bg 


eg 


57r/12 


af 


bf 


cf 


7r/3 


ae 


be 


ce 






p.|60| 


7r/4 


ad 


bd 


cd 


tt/6 


ac 


be 


cc 


7r/12 


ab 


bb 


cb 


p = 


aa 


ba 


ca 






p.^ 


= 


: 


TT 


7r 



dl 


el 


fl 


dk 


ek 


fk 


dj 


ej 


fj 


di 


ei 


fi 


p. 35 




dh 


eh 


Qr 


[dg 


eg 


fg 


df 


ef 


ff 


de 


ee 


fe 




dd 


ed 


fd 


dc 


ec 


fc 


db 


eb 


ec 


da 


ea 


fa 




P.y 


TT 


7r 




4 


3 


12 



gi 


hi 


11 


gk 


hk 


ik 


gj 


hj 


ij 


gi 


hi 


ii 


P- 


66 




gh 


hh 


ih 


[gg| 


hg 


ig 


gf 


hf 


if 


ge 


he 


ie 




p. 


62 




gd 


hd 


id 


gc 


he 


ic 


gb 


hb 


ib 


ga 


ha 


ia 




p. 
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7r 




27r 


2 


12 


3 



jl 


kl 


11 


jk 


kk 


Ik 


m 


kj 


Ij 


ji 


ki 


li 


p. 37 






kh 


Ih 




kg 


Ig 


jf 


kf 


If 


je 


ke 


le 


p.|6^ 


jd 


kd 


Id 


jc 


kc 


Ic 


jb 


kb 


lb 


ja 


ka 


la 




pN 


37r 




llTT 



algorithm described in (A.4) above. The layout of the plots is shown in the chart 
m Table |. The plots for 6 = TT or p = TT beyond the top and righ t of t his chart are 
the same as those for 6 — Q or p — Q, because of the periodicity ( A. 17 ). 

The "ultra-smooth" scaling function with m^^'''^ (z) divisible by (1 -I- 2;)'^, shown 
in Figure ^ above, and its counterpart under the symmetry (A. 18), lie at the posi- 
tions shown by a round solid point (•) in both Ta ble p| and Figure ^. 

The scaling functions from the g — 2 family in |BrJo99b appear as subsets of 
the 3 = 3 family here, supported on various subintervals of [0, 5] of length 3. The 
correspondence results, for particular values of {6, p), from the vanishing of two of 
the Oi coefficients, and the equality of the other four a^'s to the four coefficients of 
the g = 2 family. The values of {0, p) corresponding to continuous scaling functions 
in the g = 2 family | |CoHc92| , |CoHc94 |Wan95| , |Wan96| , pauL92t (see | |BrJo99bt 
Remark 3.1]) are indicated in Table IJ. The values of {9,p) that give these known 
continuous scaling functions are indicated graphically in Figure ^, along with the 
vanishing-moment points where the polynomial rOf^'''^ (z) is divisible by (1 + z)^ 
and by (1 -I- zf' (see Proposition 8.12| ), and the tight-frame cases (see Proposition 
A.l). Some regions of the (0,p) plane where the cascade approximants do not 
converge to a continuous scaling function are also indicated in the same figure. 

Putting all the 144 pictures together as illustrated in Table ||, we get graphic 
support for the observation that the two spin-vectors in the factorization ( |A.g| ) 
produce more smoothness of x 1— > Lpe,p {x) when they are not aligned, i.e., off the 
diagonal 9 = p. It also shows that, close to one of the true Haar wavelets, i.e., when 
if is the indicator function of some [k, fc -I- 1), there is a continuous ip, while close 
to a mock Haar wavelet (i.e., one that is only a tight frame) it appears that the 
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Support interval: 

X e [0, 3] 
xe [1,4] 
X e [2, 5] 



?, p) values: 



Continuous 



, (x) at: 



9 = 37r/4} 9 = 'Sir/ A, p e (0, 7r/4) U (37r/4, tt) 
p = 0} p = 0, 6le (7r/4,7r/2)U(7r/2,37r/4) 
6» = 7r/4} 6' = 7r/4, (0,7r/4)U(37r/4,7r) 




Figure 4. T/im curved lines (in the four corners): vanishing first 
moment of V' (Proposition 8.12| (a)). Round solid points (•): van- 
ishing second moment of tp; "ultra-smooth" wavelet scaling func- 
tion (Proposition |8.12| (b), Figure |2|). Thick straight lines: em- 
bedding of (continuous portion of) g — 2 family m g = 3 fam- 
ily (Table Round open points (O): translated Haar functions 
within 5 = 2 family (plots "da", "dd", "dj", "ga", "ja", "jd", 
"jj"). Square points: tight frames (Proposition A.l , plots "aa", 
"dg", "gg", "jg"). Shading: divergence due to ag > 1 or 05 > 1. 



graph of the scaling functions have Hausdorff dimension > 1, hence the "fractal" 
appearance. 

Acknowledgements. We are very grateful to Ola Bratteli, Ken Davidson, and David 
Kribs for enlightening discussions, and to Brian Treadway and Cymie Wehr for 
expert typesetting, and graphics artwork. 
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2 


2 




1 

2 


1 ^' 
1 

2 



Q123-35Q123-35Q123-35 



dd: = 7r/4, p = 7r/4 cd: 61 = 7r/3, p = 7r/4 fd: 9 = bw/U, p = w/A 



— ^ 


2 

" — ^ 


2 




1 

2 


1 

2 



Q123-35Q123-35Q123-35 



dc: 6> = 7r/4, p = 7r/6 ec: 61 = n/3, p = 7r/6 fc: 6 = 57r/12, p = 7r/6 



— 


2 


2 




1 

2 


„ V^^i^ — 

1 

2 



Q123-35Q123-35Q123-35 



db: 6l = 7r/4, p = 7r/12 eb: 61 = 7r/3, p = 7r/12 fb: 6* = 57r/12, p = 7r/12 











1 

2 


2 



Q123-35Q123-35Q123-35 



da: 61 = 7r/4, p = ea: 61 = 7r/3, p = fa: = 5tt/U, p = 
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jd: e = 37r/4, p = 7r/4 kd: 6 = 57r/6, p = 7r/4 Id: = ll7r/12, p = 7r/4 




■J 






jc: 6» = 37r/4, p = 7r/6 kc: 61 = 57r/6, p = 7r/6 Ic: 6 = llir/U, p = ir/6 


'.A 


f 


2 


1 

2 


■r 

2 


2 



jb: 6* = 37r/4, p = 7r/12 kb: 61 = 57r/6, p = 7r/12 lb: (9 = ll7r/12, p = 7r/12 



ja: 6* = 37r/4, p = ka: 61 = 57r/6, p = la: 61 = ll7r/12, p = 
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nh: e = 0, p = 77r/12 bh: = 7r/12, p = 77r/12 ch: 9 = n/6, p = 7^/12 



Q I 

af: (9 = 


2 3-35 

= 0, p = 57r/12 











'- f-r-^^- 




Q I 

bg: ^ = 


2 3 J 5 

:7r/12, p = 7r/2 


f 








Q I 

bf: 6 = 


?r/12, p'= 57r/12 ' 








■ ■ 



ae: ^ = 0, p = tt/S be: = 7r/12, p = 7r/3 ce: ^ = 7r/6, p = 7r/3 
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gh: e = 7r/2, p = 77r/12 hh: 9 = 77r/12, p = 7-k/12 ih: = 27r/3, p = I-k/12 




gg: e = 7r/2, p = 7r/2 hg: 6 = 77r/12, p = 7r/2 ig: 6> = 27r/3, p = 7r/2 




gf: e = 7r/2, p = 57r/12 hf: 61 = 77r/12, p = 57r/12 if: 9 = 27r/3, p = 57r/12 




ge: 6» = 7r/2, p = 7r/3 he: 6» = 77r/12, p = 7r/3 ie: 61 = 27r/3, p = 7r/3 
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jh: e = 37r/4, p = 7^/12 ' kh: 9 = 57r/6, p= 7^/12' Ih: 6* = ll7r/12, p = 77r/12 




jg: 6' = 37r/4, p = 7r/2 kg: 6 = 57r/6, p = 7r/2 Ig: g = ll7r/12, p = 7r/2 



■i ■ ■ 



^^1 



jf: 6* i 37r/4, p = 57rVl2 ' kf: 61 = 57r/6, p= 5n/12 ' If: 61 = ll7r/12, p = 57r/12' 



je: 6* = 37r/4, p = 7r/3 ke: 6* = 57r/6, p = 7r/3 le: 61 = ll7r/12, p = 7r/3 



64 



PALLE E. T. JORGENSEN 




al: 6» = 0, p = ll7r/12 bl: 6* = 7r/12, p = ll7r/12 cl: = ir/G, p = lln/12 





. "* . + 






ak: 6» = 0, p = 57r/6 bk: 61 = 7r/12, p = 57r/6 ck: 9 = 7r/6, p = 57r/6 











Q I 

aj: e-- 


= 0, p = 37r/4 






12 3-35 




ai: 6» = 0, p = 27r/3 bi: 6* = 7r/12, p = 27r/3 ci: 6* = tt/6, p = 27r/3 
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a -s/^^^^\/ 


2 

1 ^^^^ 




1 

2 


1 
2 



dl: 6* = 7r/4, p= ll7r/12 cl: 9 = ir/S, p = lln/12 fl: 6* = 57r/12, p = ll7r/12 





2 

1 ^^^^ 


2 

" — — 




1 

2 


2 



dk: 61 = 7r/4, p = 57r/6 ck: 61 = 7r/3, p = 57r/6 fk: 6* = 57r/12, p = 57r/6 



^^^^ 



dj: 61 = 7r/4, p = 37r/4 ej: = tt/3, p = 'Sir/A fj: = 57r/12, p = 37r/4 



31 ■ 




't- •^■ ■- -I ■ 



di: 6» = 7r/4, p=27r/3 ei: 61 = 7r/3, p = 27r/3 H: 9 = bn/U, p = 2^/3 
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gl: = -it/2, p = ll7r/12 h i: = Iw/U, p = llw/U il: = 2tt/3, p = llw/U 




gk: 6> = 7r/2, p = 57r/6 hk: = Iw/U, p = 57r/6 ik: = 2it/3, p = 57r/6 




Q123A5Q123A5Q123A5 



gj: = tt/2, p = 37r/4 hj: = Ttt/U, p = 37r/4 ij: = 2tt/3, p = 37r/4 




Q123A5Q123A5Q123A5 



gi: = Tr/2, p = 2n/3 hi: = 77r/12, p = 2tt/3 ii: = 2tt/3, p = 2tt/3 
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jl: e = 37r/4, p = IItt/u' kl: 9 = 57r/6, p = Utt/U Ih 6 =\i7t/\2, p= ll7r/12 





jk: 6» = 37r/4, p = 57r/6 ' kk: 6*'= 5n/6, p = 57r/6 ' Ik: 6* = llw/u'p = Stt/s' 




jj: 6» = 37r/4, p = 37r/4 kj: 6* = 57r/6, p = 37r/4 Ij: 6* = ll7r/12, p = 37r/4 




ji: 6* = 37r/4, p= 27r/3 ki: 61 = 57r/6, p = 27r/3 li: 61 = ll7r/12, p = 27r/3 
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